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Abstract 


The  overall  objective  of  the  project  was  to  design  and  study  new  flame  stabilizers  for  aug- 
mentor  systems  that  lead  to  improved  static  flame  stability.  For  this  purpose,  the  project  con¬ 
sisted  of  two  major  activities:  (a)  development  of  new  augmentor  flame  holder  designs  with 
improved  static  flame  stability,  and  their  experimental  characterization;  and  (b)  demonstra¬ 
tion  of  numerical  predictions  of  the  flame  behavior  using  the  resulting  experimental  designs 
and  measurements. 

The  experimental  studies  investigated  the  effect  of  cavity-based  geometries  under  the 
following  four  flame  configurations:  (i)  partially  prenrixed  flame  in  high  temperature  vitiated 
flow;  (ii)  partially  premixed  flame  at  room  temperature  air  flow;  (iii)  fully  premixed  flame; 
and  (iv)  hybrid  fully  and  partially  prenrixed  flames.  The  first  flame  condition  was  used  to 
study  the  effect  on  liftoff  heights.  Blowout  studies  were  conducted  in  the  remaining  three 
flame  conditions.  All  studies  demonstrated  improved  static  flame  stability  of  the  new  flame 
holders.  This  behavior  was  explained  using  measurements  of  species,  temperature,  and 
velocity  held. 

Based  on  the  selected  flame-holder  configurations  from  experiments,  two  computational 
modeling  studies  were  conducted  to  establish  a  new  formulation  of  Damkohler  number  sim¬ 
ilarity  and  to  assess  new  combustion  modeling  techniques  to  predict  static  stability  of  the 
flame-holder. 

Simulations  of  how  and  combustion  for  the  new  flame  holder  designs  compared  quali¬ 
tatively  well  with  the  experiments.  However,  there  were  disagreements  between  simulation 
and  experimental  results.  In  the  simulations,  the  hames  remained  attached  to  the  bluff  body 
for  all  cases  and  hames  formed  on  the  edge  of  the  bluff  body  wake.  Both  these  features 
were  not  seen  in  the  experiments.  It  was  speculated  that  the  reason  for  disagreement  was  a 
higher  level  of  inflow  turbulence  in  the  experiment  or  a  very  non-isotropic  how  induced  by 
the  turbulence  mesh.  However  further  investigation  is  required. 
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1  Introduction  and  Project  Objectives 


The  overall  objective  of  the  project  was  to  design  and  study  new  flame  stabilizers  for  augmentor 
systems  that  lead  to  an  improved  static  flame  stability  computationally  and  experimentally.  The 
main  objectives  of  the  project  were:  i)  development  of  various  newly  designed  bluff  body  (aug¬ 
mentor)  bases,  ii)  experimental  investigation  of  static  flame  stability  with  the  various  prototype 
base  geometries,  and  iii)  demonstration  of  numerical  predictions  of  the  experiments.  In  accor¬ 
dance  with  these  objectives,  the  experimental  studies  focused  on  proposing  local  cavity-based 
base  geometries  and  investigating  their  effect  on  liftoff  heights  and  blowout  limits  for  the  follow¬ 
ing  flame  configurations:  i)  a  partially  premixed  flame  in  a  high  temperature  vitiated  flow  (liftoff 
study);  ii)  a  partially  premixed  flame  in  room  temperature  air  flow  (blowout  study);  iii)  a  fully 
premixed  flame  (blowout  study);  and  iv)  hybrid  fully  and  partially  premixed  flames  (blowout 
study).  The  detailed  experimental  results  of  the  liftoff  and  blowout  studies  will  be  discussed 
in  section  2.  The  computational  studies  will  be  presented  in  section  3.  Two  different  compu¬ 
tational  approaches  were  used  to  assist  in  the  development  of  new  flame- holding  technologies. 
The  first  was  a  new  formulation  of  Damkohler  (Da)  number  similarity.  This  technique  allowed 
for  a  proper  definition  of  Da,  such  that  similarity  conditions  was  defined  for  the  assessment  of 
flame  stabilization.  With  this  method,  the  assessment  of  afterburner  static  stability  under  typ¬ 
ical  operating  conditions  could  be  performed  using  low-speed  experiments.  The  second  was  an 
assessment  of  the  combustion  modeling  techniques  to  predict  static  stability  of  the  flame-holder 
configurations  presented  in  Section  2. 

2  Flame  Holder  Development 

2.1  Experimental  Setup  and  Flame  Holder  Design 

2.1.1  Setup 

As  shown  in  Fig.  1,  the  experimental  setup  for  the  liftoff  study  consisted  of  six  parts:  i)  a  main 
premixed  burner  to  provide  post-combustion  gases  in  the  vicinity  of  a  bluff  body,  ii)  a  bluff  body 
through  which  methane  fuel  was  injected  to  form  jets  in  crossflow  with  the  high-temperature 
post-combustion  stream,  iii)  an  additional  cooling  capability  in  the  post-combustion  stream  of 
the  premixed  burner  and  upstream  of  the  bluff  body  to  control  the  vitiated  air  temperature,  iv) 
a  turbulence  grid  to  mimic  high-Reynolds  number  effects,  v)  a  stereo  PIV  system  to  measure 
the  local  free-stream  speed  and  the  velocity  fields  in  the  vicinity  of  the  bluff  body,  and  vi) 
an  intensified  camera  for  CH  chemiluminescence  detection  by  which  the  liftoff  heights  of  the 
methane  diffusion  flame  were  characterized.  In  the  following  sections,  a  detailed  description  of 
the  experimental  facility  is  provided. 
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Figure  1:  Schematic  of  the  top  view  of  the  experimental  setup 
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2.1.2  Implementation  of  Premixed  Burner 

The  premixed  burner  (on  loan  from  Professor  Robert  Dibble’s  lab  at  University  of  California, 
Berkeley)  was  located  centrally  in  a  50  cm  x  50  cm  cross-sectional  area  vertical  wind  tunnel.  The 
wind  tunnel  could  generate  3.5  m/s  flow.  The  burner  had  2184  holes,  each  having  a  diameter 
of  1/16  inch  on  the  top  surface.  The  top  surface  had  an  overall  diameter  of  8.27  inches  to  allow 
the  array  of  tiny  flames  to  generate  a  uniform  temperature  profile  over  a  wide  region.  Air  was 
delivered  by  a  1.5  hp  blower  and  was  mixed  with  methane  approximately  3  m  upstream  of  the 
burner  exit.  The  equivalence  ratio  of  the  methane/air  mixture  was  kept  fixed  at  0.64  in  the 
current  study.  Furthermore,  the  resulting  temperature  and  speed  of  the  post  combustion  gas 
were  ~1315K  and  5.7  m/s. 


2.1.3  Basic  Bluff-Body  Design 

The  key  issues  in  the  bluff  body  flame  holder  design  were  the  geometrical  shape  of  the  flame 
holder,  the  dimensions,  and  the  ability  of  the  material  to  survive  in  the  high  temperature 
environment.  The  basic  geometry  was  chosen  in  consultation  with  Dr.  Jeffrey  Lovett  from  Pratt 
&  Whitney  and  Dr.  Marios  Soteriou  from  United  Technologies  Research  Center.  As  shown  in 
Fig.  2,  the  basic  geometry  of  the  bluff  body  consisted  of  a  4.5  inches  long,  linch  thick,  and 
6  inches  wide  rectangular  body  with  a  1/2  inch  radius  rounded  nose.  A  porous  ceramic  (fused 
silica,  63%  nominal  porosity),  which  is  machinable,  has  high  service  temperature  (~2000  K), 
and  has  good  resistance  to  thermal  crack  propagation,  was  used  to  build  the  bluff  body.  The 
high  service  temperature  eliminated  the  need  for  complex  cooling.  Methane  was  delivered  from 


Figure  2:  Detailed  schematic  of  the  bluff  body. 


six  1/8  inch  (OD)  fuel  injection  nozzles  (three  on  each  side)  branched  from  a  single  3/8  inch  OD 
fuel  delivery  pipe.  The  nozzles,  which  were  located  1/2  inch  upstream  of  the  base  of  the  bluff 
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body  with  1-inch  spacing,  produced  two  sets  of  opposing  jets  in  the  crossflow.  Throughout  the 
study,  the  initial  speed  of  the  methane  jet  and  the  corresponding  square  root  of  the  momentum 
ratio  (r)  between  the  individual  methane  jets  and  the  crossflow  were  kept  fixed  at  ~3.7  m/s  and 
~0.9,  respectively.  The  origin  of  the  coordinate  system  was  placed  on  the  symmetry  planes  of 
the  bluff  body  and  1  inch  downstream  of  the  physical  base,  as  shown  in  Fig.  2.  Furthermore, 
the  dotted  volume  in  Fig.  2  would  be  filled  later  with  various  physical  inserts. 

2.1.4  Implementation  of  the  Cooling  Capability 

The  copper  water-cooled  grid,  which  has  a  geometry  similar  to  the  condenser  coil  of  a  refrigerator, 
was  placed  on  the  top  of  the  main  burner  to  control  the  temperature  of  the  vitiated  flow.  Water 
at  room  temperature  flowed  through  a  1/4  inch  (OD)  serpentine-shaped  copper  tube  that  had 
11  turns  within  a  15  inches  by  18  inches  rectangular  region.  Numerous  1/16  inch  (OD)  copper 
wires  were  welded  on  the  copper  tube  to  improve  the  heat  conduction  along  the  streamwise 
direction  of  the  water  flow.  Using  the  water  cooled  grid,  we  were  able  to  reduce  the  temperature 
of  the  post-combustion  gas  by  ~100  K  from  its  original  temperature,  which  was  at  ~1400  K. 
The  ability  to  reduce  the  temperature  was  an  important  aspect  of  this  study,  because  without 
this  reduction,  the  methane  jets  ignited  immediately  upon  exiting,  whereas  with  the  cooler 
free-stream,  the  jets  mixed  with  the  free-stream  and  only  ignited  in  the  wake  of  the  bluff  body. 

2.1.5  Implementation  of  the  Turbulence  Grid 

The  turbulence  grid  was  placed  ~1  inch  downstream  of  the  top  surface  of  the  water-cooled  grid. 
The  water-cooled  grid  was  composed  of  two,  top  and  bottom,  layers  of  1/2  inch  OD,  8  inch 
long  ceramic  (99.9  %  AI2O3)  tube  arrays.  Each  array  was  formed  by  6  parallel-placed  ceramic 
tubes  with  1.25  inches  inter-distance  (36%  open  area).  The  top  and  bottom  arrays  were  placed 
perpendicular  to  each  other  such  that  the  open  area  became  an  exact  square.  Using  the  grid, 
we  could  increase  the  freestream  turbulent  fluctuation  from  ~5  %  (with  no  grid)  to  8  %  (with 
the  grid). 

2.1.6  Implementation  of  Stereo  PIV 

For  the  PIV  studies,  3  //m  nominal  diameter  AI2O3  particles  were  seeded  into  the  flow,  and  a 
«100  mJ/pulse  (532  nm)  double-pulsed  Nd:YAG  laser  (New  Wave,  Gemini  PIV)  with  «1  mm 
sheet-thickness  was  used  to  illuminate  the  particles.  The  AI2O3  particles  were  mixed  densely 
with  the  methane  and  air  streams  such  that  adequate  seeding  density  could  be  achieved  even  in 
the  high  temperature  environment.  The  resulting  Mie  scattering  was  detected  by  two  double¬ 
exposure  CCD  cameras  (La  Vision,  Flow  Master)  mounted  with  a  70°  included  angle. 


2.1.7  CH  Chemiluminescence  Detection 


The  mean  and  fluctuating  components  of  the  flame  liftoff  height  were  used  as  a  criterion  for  the 
performance  of  various  bluff  body  geometries  tested  in  the  current  study.  Intensified  and  filtered 
imaging  of  CH  chemiluminescence  was  carried  out  to  measure  the  liftoff  height.  For  this  imaging, 
an  intensified  camera  (Princeton  Instrument,  PI-MAX)  with  a  gate  width  of  100  /xs  and  a  ~4305 
nm  band  pass  filter  was  used.  Representative  examples  of  the  instantaneous  chemiluminescence 
images  are  provided  in  Fig.  3.  The  white  dotted  lines  represent  the  instantaneous  location  of  the 
flame  base.  From  each  image,  the  flamebase  location  was  determined  manually  and  the  obtained 
locations  from  200  images  were  ensemble  averaged. 


56.5  mm 

< - ► 


Figure  3:  Representative  instantaneous  images  of  intensified  chemiluminescence.  Gate  width 
was  100  ms.  White  dotted  lines  represent  the  instant  locations  of  the  flamebase. 


2.2  Lift-Off  Experiment 

2.2.1  Preliminary  Geometric  Modifications  of  the  Bluff  Body  Base 

The  base  geometry  of  the  bluff  body  was  modified  by  addition  of  blocks  with  different  shapes 
to  the  dotted  region  in  Fig.  2.  Nine  7  inch  long  blocks  with  different  cross-sections,  as  shown 
in  Fig.  4,  were  used.  The  blocks  were  categorized  with  respect  to  their  shape  as  Reference  (R), 
V  cavity,  D  cavity,  S  or  R+S  cavity,  and  T  and  I  blocks.  The  Reference  case  (R)  was  a  solid 
block  extension  with  one  inch  square  cross  section.  The  R  block  modified  the  geometry  by 
increasing  the  height  of  the  bluff  body  to  5.5  inches.  The  fuel  injection  holes  were  unchanged, 
and,  therefore,  the  Reference  case  modified  the  setup  by  moving  the  fuel  injection  jets  in  the 
upstream  direction  by  1  inch.  The  Reference  case  was  used  to  form  the  V  cavity  block  by 
cutting  a  0.43  inch  deep  V  shape  with  angle  60°  along  the  block  (see  Fig.  4b  for  its  detailed 
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Figure  4:  Schematics  of  various  blocks,  (a)  Reference  (or  R),  (b)  V  cavity,  (c)  D  cavity,  (d)  15°  S 
cavity,  (e)  25°  S  cavity,  (f)  90°  S  cavity  (or  R+S),  (g)  T,  (h)  I,  and  (i)  15°  S  cavity-out-of-phase 
geometries. 
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dimensions).  Similarly,  the  D  cavity  block  was  constructed  by  cutting  a  0.3  inch  deep,  0.5  inch 
wide  rectangular  shape,  as  shown  in  Fig.  4c.  The  S  cavities  had  a  more  complex  shape,  formed 
from  the  Reference  case  by  cutting  0.5  inch  wide  individual  side-wise  cavities  with  15°  (Figs.  4d 
and  4i),  25°  (Fig.  4e),  and  90°  (Fig.  4f)  slopes,  respectively.  (The  inclined  design  was  motivated 
by  the  fastback  designs  used  in  automobile  designs,  where  the  flow  pattern  changed  based  upon 
inclination  angle.  For  details  see  Sovran  et  al.  (1978)).  The  cavities  were  1  inch  apart  and  were 
in  phase  with  the  three  fuel  nozzles  for  cases  d,  e,  and  f,  whereas  case  i  was  designed  to  be  out  of 
phase.  The  effect  of  in-phase  and  out-of-phase  geometries  will  be  discussed  later.  Furthermore, 
the  T  block  (Fig.  4g)  had  an  equilateral  triangular  shape  with  1-inch  height  and  1-inch  base 
length.  Addition  of  the  T  block  resulted  in  two  symmetric  recirculation  zones.  Finally,  as  shown 
in  Fig.  4h,  the  I  geometry  had  a  0.25  inch  x  1  inch  rectangular  cross  section. 


2.2.2  Measurement  of  Flame  Liftoff  Heights 


Figure  5  shows  the  ensemble-averaged  flame  liftoff  height  for  each  different  case  discussed  in 
section  2.2.1.  For  all  cases  the  temperature  of  the  post-combustion  gas,  the  cross  flow  speed, 


Figure  5:  Liftoff  heights  for  various  bluff  body  base  designs.  For  all  cases,  the  temperature 
and  the  speed  of  the  crossflow  and  the  square  root  of  the  momentum  ratio  between  the  high- 
temperature  post-combustion  gas  and  the  injected  methane  are  kept  identical  as  1300  K,  5.7 
m/s,  and  «0.9,  respectively. 


and  the  square  root  of  the  momentum  ratio  between  high-temperature  post-combustion  gas  and 
injected  methane  were  kept  identical  at  1300  K,  ~5  m/s,  and  ~1,  respectively.  The  equivalence 
ratio  of  the  main  burner  was  0.63.  The  liftoff  heights  were  determined  by  averaging  individual 
liftoff  heights  from  200  randomly  sampled  CH  chemiluminescence  images  of  the  flame.  The 
liftoff  height  of  the  V  cavity  (36.2  mm)  showed  very  marginal  improvement  in  comparison  with 
that  of  the  Reference  case  (~36.3  mm),  while  the  D  cavity  (33.9  mm)  showed  a  more  noticeable 
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decrease  of  the  liftoff  height.  The  most  significant  improvement,  however,  was  seen  in  the  cases 
of  the  three  S  cavities  (22.4  mm,  22.1  mm,  and  22.5  mm  for  15°,  25°,  and  90°  wall  slopes, 
respectively)  and  the  T  (22.8  mm)  and  I  (22.3  mm)  geometries.  The  flame  liftoff  heights  were 
not  a  strong  function  of  the  slope  of  the  S  cavities  in  the  range  investigated  here  even  though 
their  recirculation  patterns  and  vortex  structures  could  be  noticeably  different  from  each  other 
(Sovran  et  al. ,  1978). 

2.2.3  Design  Optimization  (D+S,  T+S  Cavities) 

Based  on  the  results  from  the  geometries  discussed  in  Section  2.2.1,  we  refined  our  designs  to 
benefit  from  the  presence  of  both  the  D  cavity  (re-entrant  shape)  and  the  90°  S  cavity,  which 
introduced  local  three-dimensionality  to  form  a  so  called  D+S  cavity  block  shown  in  Fig.  6a. 
The  D  cavity  part  was  ~0.5  inch  deep  and  ~0.5  inch  wide,  and  the  90°  S  cavity  is  0.5  inch  wide, 
1  inch  high,  and  0.4  inch  deep.  The  T  and  the  90°  S  cavity  were  combined  to  build  a  so  called 
T+S  block  shown  in  Fig.  6b.  The  S  cavity  was  0.5  inch  wide  and  0.25  inch  deep  (Shallow  Cut 
T+S).  The  flame  liftoff  height  measured  in  the  identical  flow  condition  as  in  Section  2.2.2  and 


Figure  6:  Schematics  of  (a)  D+S  cavity  and  (b)  T+S  geometries. 


the  cavity  placement  (in  phase  with  the  methane  nozzle)  with  that  of  Fig.  5  was  17.4  mm  in  both 
cases,  a  52  %  liftoff  height  decreased  in  comparison  with  the  Reference  case.  The  flamebase  in 
these  cases  mostly  resided  in  the  vicinity  of  (sometimes  inside  of)  the  S  cavities,  which  resulted 
in  less  fluctuation  of  the  flamebase  (the  standard  deviations  of  the  flamebase  were  ~7.5  mm  for 
D+S  block  and  7  mm  for  T+S  block,  respectively)  when  compared  to  that  of  the  Reference  case 
(~13  nun).  We  believed  that  local  back  steps  (see  Fig.  6)  formed  by  the  intersection  of  the  D 
and  S  cavities  and  the  T  and  S  cavities  might  provide  a  desirable  velocity  field  to  maximize  the 
flame  stability,  as  well  as  the  presence  of  the  S  cavity,  which  confined  the  injected  methane  to 
a  more  localized  cavity  region  and  prevents  the  methane  jets  from  becoming  excessively  diluted 
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by  the  surrounding  post-combustion  gas.  Such  a  locally  fuel-rich  (compared  to  its  surrounding) 
cavity  might  act  like  a  piloting  region  that  assists  in  better  holding  of  the  overall  downstream 
flame.  It  was,  therefore,  likely  that  this  cavity  fuel  concentration  effect  would  not  be  directly 
transferable  to  a  flame  configuration  such  as  a  fully  premixed  flame  with  a  streamlined  bluff 
body  under  far-upstream  fuel  injection  that  had  been  studied  more  recently  in  the  literature 
(Shanbhogue  et  al,  2009&).  A  more  detailed  description  of  the  fuel  concentration  field  will  be 
discussed  in  the  blowout  study  portion  of  the  current  report. 

2.2.4  Dependence  of  the  Liftoff  Heights  on  Cavity  Depths  and  Placements 

Figure  7a  shows  the  liftoff  height  variation  of  the  D+S  cavity  block  using  three  different  D 
cavity  depths  (0  inch,  0.5  inch,  and  0.8  inch),  whereas  Fig.  7b  shows  the  effect  of  different  cavity 
placements  with  respect  to  the  methane  nozzle  location,  i.e.,  in-phase  or  out-of-phase  cavities. 
Furthermore,  Fig.  7  demonstrates  the  liftoff  height  dependence  on  the  S  cavity  depth  for  the 
T+S  geometry.  In  all  three  cases,  the  flow  condition  was  identical  to  that  of  Section  2.2.2.  From 
Fig.  7a  the  liftoff  height  for  a  moderate  D  cavity  depth  (17.4  mm)  was  significantly  lower  than 
that  of  the  other  two  cases  (22.5  mm  in  the  0  inch  case  and  33.2  mm  in  the  0.8  inch  case). 
Neither  the  significantly  shallow  (0  inch  case)  nor  the  significantly  deep  D  cavity  depth  (0.8  inch 
case)  were  optimal  for  flame  stability.  Hence,  a  D  cavity  with  unity  aspect  ratio  was  better. 

The  effect  of  phase  matching  between  the  S  cavity  and  the  methane  injection  nozzle  also 
played  a  central  role  in  liftoff  height,  as  shown  in  Fig.  7b.  The  liftoff  height  increased  by  more 
than  90  %  by  shifting  the  D+S  cavity  by  a  half  period  such  that  methane  fuel  was  injected  into  a 
no-cavity  region.  The  importance  of  the  phase  matching  was  also  clear  when  different  geometries 
were  used.  For  example,  the  liftoff  height  of  the  completely  in-phase  15°  S  cavity  (Fig.  4d)  was 
~22.5%  lower  than  that  of  the  half-in-phase/half-out-of-phase  15°  S  cavity  (Fig.  4i).  Addition¬ 
ally,  the  flame  liftoff  height  was  a  function  of  the  depth  of  the  S  cavity,  as  shown  in  Fig  7c.  The 
liftoff  height  from  a  moderate  S  cavity  depth  (0.25  inch)  showed  a  lower  value  than  that  of  the 
0  and  0.375-inch  depth  cases. 

2.2.5  Summarizing  Plot  of  Liftoff  Height  Measurements 

As  shown  in  the  summarizing  plot  (Fig.  8),  the  geometries  containing  local  cavities  (S  cavities) 
with  moderate  aspect  ratio  and  depth,  especially  the  D+S  and  the  T+S  cavities  placed  in-phase 
with  methane  injection,  showed  over  a  50  %  decrease  in  the  liftoff  heights  in  comparison  with  the 
Reference  case.  In  summary,  the  results  shown  here  suggested  that:  (a)  a  combination  of  two- 
and  three-dimensional  changes  to  the  base  region  (D+S  and  T+S)  provided  the  lowest  liftoff 
height;  and  (b)  the  alignment  of  the  local  cavities  with  the  jet  in  crossflow  resulted  in  optimal 
performance. 
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Figure  7:  Liftoff  height  dependence  on  (a)  the  cavity  depth  (D  cavity  depth  of  the  D+S  cavity), 
(b)  the  placement  of  the  bluff  body  base  (D+S  case),  and  (c)  the  S  cavity  depth  (T+S  case). 


(c) 


0"  depth  0.25"  depth  0.375"  depth 


14 


Reference  D+S  Cavity  (Aspect  T+S  (1/4"  Shallow 
Ratio  =  -1)  Cut) 


Figure  8:  Minimum  liftoff  heights  observed  in  the  current  study.  The  D+S  geometry  with 
moderate  (~1)  aspect  ratio  and  the  T+S  geometry  with  0.25  inch  S  cavity  depth  showed  the 
minimum  liftoff  heights  (~50%  of  the  Reference  case). 


2.2.6  Velocity  Fields  Obtained  from  Stereo  PIV 

The  six  planes  investigated  in  the  stereo  PIV  measurements  are  illustrated  in  Fig.  9.  Three 
planes  (Pxy_C,  Pxy_Q  and  Pxy_H)  were  parallel  to  xy  plane,  where  x  was  a  coordinate  parallel 
to  the  initial  direction  of  methane  jet  injection,  and  y  represented  a  spanwise  coordinate  (see 
Fig.  2).  The  other  three  (Pyz_C,  Pyz_Q,  and  Pyz_H)  were  parallel  to  yz  plane.  The  Pxy_C 
plane  passed  through  the  center  of  spanwise  (S)  cavities  in  the  D+S  or  T+S  geometries,  and  the 
Pxy_Q  and  Pxy_H  planes  were  offset  0.25  inch  and  0.5  inch  from  the  center  plane,  respectively. 
Hence,  the  plane  Pxy_Q  passed  through  the  side  wall  of  the  S  cavity,  and  Pxy  H  was  located 
at  the  center  of  the  no-cavity  region.  The  planes  Pyz_C,  Pyz_Q,  and  Pyz  H  were  located  at  the 
center  of  the  bluff  body,  0.25  inch  from  the  center  and  0.5  inch  from  the  center,  respectively.  The 
Pyz_H  plane  intersected  the  edge  of  the  bluff  body.  Throughout  this  section,  only  the  reference 
geometry,  the  D+S,  and  the  T+S  geometries  with  the  lowest  liftoff  heights  (see  Fig.  8)  were 
considered. 

Figures  10a,  b,  and  c  show  ensemble-averaged  PIV  images  from  100  instantaneous  snapshots 
for  the  Reference,  D+S,  and  T+S  geometries  along  the  Pxy_C  plane,  respectively.  Figure  lOd 
shows  the  corresponding  centerline  velocities  along  the  y  coordinate.  The  color  coding  in  all 
the  PIV  images  provided  in  this  study  represents  the  out-of-plane  velocity  component.  For 
example,  the  red,  green,  and  blue  colors  represent  out-of  (+z  direction),  within  (in  plane  xy), 
and  into-the-paper  (-z  direction)  components  of  the  velocities,  respectively.  The  white  lines  are 
pseudo-streamlines  obtained  from  in-plane  velocity  components. 

As  shown  in  Figs.  10a,  10b,  and  10c,  xy-plane  recirculation  regions  were  present  in  all 
three  cases,  as  expected;  however,  noticeable  differences  in  terms  of  location  of  vortices  and  the 
strength  of  the  recirculation  in  the  vicinity  of  the  bluff  body  were  observed  among  the  three 
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Figure  9:  Schematic  of  the  six  planes  of  interest  in  the  PIV  measurements. 


cases.  The  center  of  vortices  (denoted  by  uzl  in  this  study)  were  further  away  from  the  bluff 
body  in  the  case  of  the  Reference  (~23  mm  from  the  base)  compared  to  that  of  the  D+S  (~14 
mm)  and  the  T+S  (~5  mm)  cases.  Furthermore,  the  minimum  vertical  velocity  in  Fig.  lOd  was 
significantly  higher  for  the  Reference  geometry  (dashed  line,  ~0.1  m/s)  than  that  of  the  D+S 
(dotted  line,  ?s-0.6  m/s)  and  the  T+S  (solid  line,  «-0.8  m/s)  geometries.  These  facts  showed 
that  the  small  distance  of  the  vortices  from  the  bluff  body  and  the  intense  near-base  recirculation 
in  the  D+S  and  T+S  geometries  caused  the  reductions  in  liftoff  heights. 

The  recirculating  regions  of  all  three  cases  generally  shared  common  features  for  the  out-of¬ 
plane  velocity  component  (w)  in  the  xy  plane,  which  are  represented  by  the  yellow  or  red  color 
in  the  recirculation  regions  of  Figs.  10a,  10b,  and  10c.  These  biased  out-of-plane  velocities  most 
likely  were  generated  due  to  the  imperfect  symmetry  in  the  current  experiment.  The  postulated 
reason  for  the  out-of-plane  velocities  was  even  more  plausible  when  one  considered  the  study  by 
Smith  &  Mungal  (1998)  in  which  they  showed  that  a  single  jet  in  crossflow  leads  to  asymmetric 
development  downstream.  Thus,  it  would  be  difficult  to  have  six  such  jets  interacting  with  the 
bluff  body  base  and  produced  a  perfectly  symmetric  flow  (the  asymmetry  was  observed  more 
clearly  in  the  yz  plane  PIV  measurements,  which  will  be  discussed  below).  A  more  careful 
inspection  of  the  results  in  Fig.  10  revealed  that  the  w  component  in  the  recirculation  region 
of  the  D+S  cavity  showed  noticeable  gradients  along  both  the  y  coordinate  (from  the  pale  blue 
at  the  bottom  to  the  red  at  y  ~16  mm)  and  the  x  coordinate  (from  the  red  at  the  center  to 
the  blue  at  x  ~  ±5  mm)  compared  to  the  Reference  case  (constant  yellow).  These  gradients 
implied  that  there  existed  out-of-plane  vortical  motions  (denoted  by  ux  1  and  ujy  1  in  directions  y 
and  x,  respectively  and  depicted  in  Fig.  10b  by  thick  white  lines)  in  the  D+S  cavity,  which  was 
one  of  the  key  elements  for  decreasing  liftoff  height  reported  in  this  study.  Those  vortices  might 
also  be  present  for  the  T+S  case,  but  they  were  not  detected.  The  plane  Pxy_C  was  located 
significantly  further  from  the  top  of  the  S  cavity  portion  in  the  T+S  geometry  than  in  the  D+S 
geometry,  which  resulted  in  lower  strength  of  the  significant  vortical  motions. 
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Figure  10:  Ensemble  averaged  PIV  images  of  (a)  the  Reference,  (b)  the  D+S,  and  the  (c)  T+S 
geometries  in  the  Pxy_C,  and  (d)  corresponding  centerline  velocities  along  the  y  coordinate.  The 
color  coding  represents  out-of-plane  velocity  component,  while  the  white  lines  show  the  pseudo 
streamlines  and  the  sketches  of  vortices. 
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Figures  11a  and  lib  show  the  ensemble- averaged  PIV  images  in  the  plane  Pxy_Q  of  the  D+S 
and  the  T+S  geometries,  respectively.  Experimental  conditions  were  identical  to  those  of  Fig.  10 
except  that  the  plane  of  interest  was  offset  by  0.25  inch  such  that  the  laser  sheet  illuminated  the 
boundary  of  the  cavity  and  no-cavity  regions.  Compared  to  the  reference  center  plane  (Pxy_C), 
the  flow  patterns  were  significantly  different  qualitatively  in  the  current  plane  for  the  D+S  cavity 
than  those  observed  for  the  T+S  cavity.  For  the  D+S  cavity,  the  strong  recirculating  flow  along 
the  negative  y-direction  was  diminished,  and  more  intense  x-directional  flow  was  observed.  This 
x-directional  motion  also  was  confirmed  by  a  subsequent  tuft  flow  visualization  in  a  non-reacting 
environment,  where  a  single  tuft  was  placed  at  the  end  of  a  long  slender  rod  to  probe  the  base 
flow.  In  accordance  with  the  tuft  visualization  experiment,  there  existed  significant  flow  motion 
from  one  S  cavity  to  the  adjacent  (or  oppositely  located)  S  cavities  through  the  D  cavity  inside. 
The  x-directional  flow  detected  here,  therefore,  was  a  consequence  of  the  cavity-cavity  flow.  For 
the  w  velocity  component  in  this  geometry,  more  intense  (bright  red)  outward  (+z  direction) 
flow  was  observed  in  the  current  plane  compared  with  that  of  the  center  plane.  This  high  flow 
intensity  was  caused  by  i)  the  current  plane  located  at  the  boundary  of  the  cavity  and  the  no¬ 
cavity  regions  was  intrinsically  asymmetric  and  ii)  there  existed  additional  inter-cavity  flows, 
from  one  cavity  to  the  adjacent  cavities  on  the  same  side.  The  more  intense  outward  flow  in  this 
plane  also  was  observed  in  the  T+S  geometry,  as  shown  in  Fig.  lib;  however,  no  other  noticeable 
difference  in  the  flow  pattern  was  observed  for  the  T+S  geometry  in  the  current  plane  compared 
with  that  in  the  center  plane. 
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Figure  11:  Ensemble  averaged  PIV  images  of  (a)  the  D+S  and  (b)  the  T+S  geometries  in  the 
Pxy-Q. 


Figures  12a  and  12b  show  the  PIV  images  at  the  center  of  the  no-cavity  regions  (Plane 
Pxy_H)  of  the  D+S  and  the  T+S  geometries,  respectively,  under  flow  conditions  identical  to 
those  of  Fig.  10.  In  the  D+S  cavity  case,  the  width  of  the  recirculation  region  of  the  D+S 
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cavity  in  the  current  plane  was  noticeably  larger  compared  to  that  in  the  plane  Pxy_C  shown  in 
Fig.  10b.  The  larger  recirculation  was  observed  because  the  stream  in  the  current  plane  did  not 
experience  the  slender  S  cavity  portion  but  instead  flowed  through  the  wider  no-cavity  region. 
All  other  flow  features  in  this  plane  for  the  D+S  case  were  similar  to  those  observed  in  Fig.  10b, 
such  as  i)  the  intense  near-base  recirculation,  ii)  the  presence  of  the  ojx  1  and  uyl  (observed  from 
the  gradient  of  color  coding  along  the  y  and  x  coordinates,  respectively),  as  well  as  w2l.  A 
more  detailed  inspection,  however,  showed  that  there  existed  quantitative  differences  between 
the  flow  fields  of  the  center  plane  (plane  Pxy_C)  and  the  current  plane,  e.g.,  the  location  of  w-1 
(~20  mm)  in  the  current  plane  was  higher  (~15  mm).  Also,  the  centerline  velocity  along  the  y 
coordinate  of  the  current  plane  was  ~20  %  smaller  than  that  observed  in  the  plane  Pxy_C. 

These  differences  also  were  observed  in  the  case  of  the  T+S  geometry  (Fig.  12b).  For  the 
T+S  geometry,  the  location  of  loz  1  was  higher  (~10  -  13  mm)  than  that  of  the  center  plane 
(~5  mm),  and  the  centerline  velocity  in  the  current  plane  was  reduced  by  30%  (not  shown)  in 
comparison  with  that  of  the  center  plane.  We  believed  that  the  higher  location  of  ajzl  and  the 
reduced  strength  of  the  recirculating  flow  were  partially  responsible  for  the  higher  liftoff  height 
in  the  case  of  out-of-phase  nozzle  placement  corresponding  to  the  nozzle  geometry  described  in 
Fig.  7b. 


Figure  12:  Ensemble  averaged  PIV  images  of  (a)  the  D+S  and  (b)  the  T+S  geometries  in  the 
Pxy_H. 


Figure  13  shows  the  ensemble-averaged  PIV  images  of  the  Reference  case  (Fig.  13a),  the 
D+S  cavity  (Fig.  13b),  and  the  T+S  (Fig.  13c)  geometry  for  the  plane  Pyz_H,  which  was 
offset  from  the  center  plane  (yz  plane)  by  1/2  inch.  The  red  and  blue  colors  represent  the 
outward  fluid  motion  (-x  direction  for  the  current  plane)  and  inward  motion  (+x  direction), 
respectively.  As  expected,  most  of  the  in-plane  (yz  plane)  velocity  vectors  uniformly  pointed 
downstream  in  all  three  cases  since  this  plane  was  very  close  to  the  freestream.  With  a  careful 
investigation  of  the  color  coding,  one  could  find  the  presence  of  the  uzl  (observed  in  Figs.  10 
and  12)  for  the  Reference  and  the  D+S  geometries  from  their  upstream  yellow  color  and  their 
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blue  color  downstream.  Another  observation  from  the  color  coding  was  that  the  outward  flow 
in  the  vicinity  of  the  no-cavity,  and  also  inward  flow  in  the  vicinity  of  the  cavity,  of  the  D+S 
geometry  was  noticeable,  unlike  the  Reference  geometry.  This  inward/outward  flow  was  because 
the  cavity  region  of  the  D+S  cavity  geometry  was  retarded  by  0.4  inch  from  the  no-cavity  region. 
The  presence  of  the  out-of-plane  velocity  also  was  found  in  the  T+S  geometry,  which  could  be 
confirmed  from  the  color  coding  gradients  in  Fig.  13c.  More  significant  inward  flow  (deep  blue 
or  black)  was  observed  in  the  no-cavity  region  downstream  of  the  T+S  in  comparison  to  that  in 
the  cavity  region,  which  was  opposite  to  the  D+S  case.  This  significant  inward  flow  was  because 
the  no-cavity  surface  of  the  T+S  is  25°  slanted  backward  from  the  current  PIV  plane,  while  the 
back  surface  of  the  S  cavity  was  parallel  to  the  plane,  which  mitigated  an  intense  inward  flow 
at  the  cavity  downstream. 

The  ensemble- averaged  PIV  images  in  the  plane  Pyz_Q  are  shown  in  Fig.  14.  Unlike  for  the 
previous  plane,  the  overall  flow  pattern  of  this  plane  did  not  point  in  the  +y  direction,  but  was 
biased  to  the  right.  This  asymmetry  probably  was  caused  by  i)  a  slight  misalignment  between 
the  burner  and  the  bluff  body,  ii)  the  complicated  interaction  between  the  crossflow  and  the 
methane  jet  (Smith  and  Mungal  1998),  and  iii)  possible  asymmetric  suction  flow  along  the  z 
direction  at  the  edges  of  the  bluff  body  due  to  the  influence  of  the  shear  layers  formed  by  the 
speed  difference  between  the  post  combustion  gas  (~5  m/s)  and  the  surrounding  non-reacting 
air  stream  (~3.5  m/s).  This  observation  was  consistent  with  the  red  shift  observed  in  Fig.  10. 
The  typical  speeds  of  the  biased  flow  measured  in  this  plane  were  quite  small,  ~0.03  m/s,  for 
the  Reference  case  and  somewhat  larger,  ~0.3  m/s,  for  the  cavity  cases. 

For  the  Reference  geometry  (Fig.  14a),  the  presence  of  the  u+1  observed  in  the  previous  planes 
was  seen  more  clearly  from  the  upstream  yellow  color  coding  with  the  subsequent  downstream 
blue.  Also,  the  low  downward  flow  speed  (flow  along  -y  direction)  confirmed  the  previous 
observations  shown  in  Fig.  10a.  For  the  D+S  geometry  shown  in  Fig.  14b,  a  more  complex 
flow  field  due  to  the  inward  (yellow)  and  outward  (blue)  flows  in  the  vicinity  of  the  cavity  and 
no-cavity  regions,  respectively,  was  seen  more  clearly.  Also,  this  geometry  had  vortices  along 
all  three  coordinates.  Along  the  x  coordinate,  the  u+1  was  present  at  the  left  bottom  corner 
of  the  figure,  which  resulted  in  a  strong  recirculation  into  the  cavity  region.  Furthermore,  one 
could  find  that  another  set  of  in-plane  vortices  (so  called,  wx2),  located  slightly  lower  than  u+1 
(see  the  figure  for  its  location),  was  present,  probably  induced  by  a  back  step  formed  by  the 
side  surface  of  the  S  cavity  and  the  bottom  surface  of  the  D  cavity.  Along  the  z  coordinate, 
the  t+1  was  present  as  could  be  conjectured  from  the  yellow  no-cavity  region  upstream  and  the 
blue  region  downstream.  Along  the  y  coordinate  the  presence  of  the  uyl  swirling  motion  from 
a  cavity  to  the  left  corner  of  the  adjacent  cavity  on  its  right  could  be  found  also  from  a  careful 
examination  of  the  color  coding  (see  the  sketch  in  Fig.  14b).  In  general,  we  believed  that  these 
multiple,  three  dimensional  small  scale,  coherent  vortices  due  to  the  presence  of  the  D  and  the 
S  cavities  mutually  dissipated  one  another.  As  a  result,  these  vortices  played  a  central  role 
in  decreasing  the  liftoff  height,  especially  in  comparison  with  the  Reference  case,  which  only 
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Figure  13:  Ensemble  averaged  PIV  images  of  (a)  the  Reference,  (b)  the  D+S  and  (c)  the  T+S 
geometries  in  the  Pyz_H. 
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Figure  14:  Ensemble  averaged  PIV  images  of  (a)  the  Reference,  (b)  the  D+S  and  (c)  the  T+S 
geometries  in  the  Pyz_Q. 
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had  the  coherent  vortex,  uzl.  These  three-dimensional  vortical  structures  might  be  expected  in 
the  T+S  geometry  (Fig.  14c),  but  they  were  not  seen  clearly  in  the  current  plane,  which  was 
probably  due  to  the  longer  distance  between  the  S  cavity  and  the  plane  of  interest. 

In  Fig.  15,  ensemble-averaged  PIV  images  in  the  yz  center  plane  (Pyz_C)  of  the  Reference 
(Fig.  15a),  the  D+S  (Fig.  15b),  and  the  T+S  (Fig.  15c)  geometries  are  shown.  Again,  the  fluid 
motions  in  the  recirculation  zone  (y  ~7  -  30  mm  for  the  Reference,  ~7  -  22  mm  for  the  D+S, 
and  ~7  -  15  mm  for  the  T+S)  of  all  three  cases  generally  pointed  to  the  right  (+z  direction)  of 
the  image,  which  was  consistent  with  Fig.  14.  Also,  other  flow  patterns  observed  in  the  previous 
planes  could  be  confirmed  in  this  plane.  For  example,  the  Reference  had  a  relatively  organized 
color  coding  pattern,  a  yellow  colored  band  7-17  mm  and  a  subsequent  blue  colored  band  more 
downstream,  which  corresponded  to  the  cuzl  vorticity.  The  complex  flow  pattern  in  the  D+S 
cavity  was  shown  too,  e.g.,  i)  outward  (yellow)  and  inward  (blue)  flows  observed  in  the  vicinity  of 
the  no-cavity  and  the  cavity  regions,  respectively,  ii)  the  presence  of  in-plane  (yz  plane)  vortices 
in  the  vicinity  of  the  boundary  of  the  cavity /no-cavity  regions  as  a  consequence  of  ujx  1,  and,  as 
a  result,  iii)  intense  downward  flows  induced  from  the  no-cavity  to  the  cavity  regions,  especially 
at  a  height  of  7  -  15  mm,  which  was  represented  as  the  intense  y-directional  downward  flow 
in  Fig.  10.  Again,  it  appeared  that  the  intense  recirculating  flows  of  this  geometry  assisted 
in  the  decrease  of  the  liftoff  height  by  promoting  localized,  but  adequate,  mixing  between  the 
cavity  flow  (primarily  fuel)  and  no-cavity  flow  (primarily  air),  which  was  beneficial  in  holding 
the  downstream  lifted  jet  flame.  In  the  case  of  the  T+S  geometry,  one  still  could  find  the 
downwards  in-plane  flows  to  the  cavity  region  (7  -  12  mm  height)  probably  due  to  the  ojx1  (not 
shown  in  the  current  region  of  interest);  however,  the  out-of-plane  (x  coordinate)  flow  motion 
of  this  geometry,  seen  from  the  color  coding,  seemed  somewhat  less-correlated  with  the  cavity 
locations  than  that  observed  in  the  D+S  cavity.  Again,  the  less  correlated  out-of-plane  flow 
pattern  seemed  to  be  caused  by  the  longer  inter-distance  between  the  S  cavity  portion  of  the 
T+S  geometry  and  the  current  plane  Pxy_C,  and,  as  a  result,  the  effect  of  the  S  cavity  was 
diminished. 

2.2.7  Vorticity  and  Strain  Rate  Fields 

The  mean  vorticity  and  mean  strain  rates  are  shown  in  Figs.  16  and  17,  respectively.  These 
fields  were  presented  for  three  cases,  Reference,  D+S  and  T+S,  in  the  xy  (Figs.  16a  and  17a) 
and  in  the  yz  planes  (Figs.  16b  and  17b).  Both  the  strain  and  vorticity  fields  exhibited  similar 
characteristics  in  the  planes  reported.  In  the  xy  plane,  the  strain  and  vorticity  were  more 
diffused  close  to  the  body  for  the  modified  geometries.  In  the  yz  plane,  the  high  shear  and 
vortical  regions  of  the  flow  were  observed  to  emanate  from  the  spanwise  geometrical  variations 
and  most  likely  played  a  role  in  the  mean  flame  behavior.  The  vortical  patterns  will  be  discussed 
further  in  the  next  section. 
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Figure  15:  Ensemble  averaged  PIV  images  of  (a)  Reference,  (b)  D+S  and  (c)  T+S  geometries 
in  the  Pyz_C. 
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Figure  16:  Mean  vorticity  fields  of  Reference,  D+S,  and  T+S  geometries  in  planes  (a)  Pxy_C 
(b)  Pyz_Q. 
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Figure  17:  Mean  strain  rate  fields  of  Reference,  D+S,  and  T+S  geometries  in  planes  (a)  Pxy_C 
(b)  Pyz_Q. 
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2.2.8  Summarizing  Schematics  of  Flow  Patterns  with  Local  Cavity  Geometries 

This  section  summarizes  the  flow  patterns  observed  thus  far.  The  flow  field  schematics  in  the 
vicinity  of  the  three  geometries  are  provided  in  Fig.  18.  The  flow  schematic  in  the  bluff  body 
downstream  of  Fig.  18  was  constructed  using  the  PIV  results  presented  in  the  previous  section, 
whereas  the  flow  inside  the  cavities  was  conjectured  based  on  a  subsequent  visualization  of  a 
tuft  flow  experiment  in  a  non- reacting  environment. 

As  shown  in  Fig.  18a,  there  exists  only  one  vortex  pair,  loz  1,  in  the  Reference  geometry. 
For  an  ideally  symmetric  experiment,  this  large  scale  vortex  would  share  identical  shedding 
frequency  with  adjacent  vortex  pairs  along  the  z  coordinate,  which  did  not  seem  beneficial  for 
improving  flame  stability  (Kurimoto  et  al.  (2004);  Paschereit  et  al.  (2006)) 

In  contrast,  multi-dimensional  vortical  structures  were  formed  along  the  z,  x,  and  y  coor¬ 
dinates  in  both  D+S  and  T+S  geometries  (Figs.  18b  and  18c,  respectively).  In  the  case  of 
the  D+S,  aside  from  uz  1,  two  additional  z-coordinate  vortical  structures  were  observed:  (a) 
cj22,  induced  by  a  back  step  formed  between  the  S-cavity  back  surface  and  the  D-cavity  bottom 
surface,  and  (b)  ujz 3,  induced  by  the  back  step  formed  by  the  bluff  body  side  surface  and  the 
S-cavity  bottom  surface.  These  vortices  had  different  shedding  locations  and  frequencies  due 
to  the  difference  in  effective  height  and  thickness  of  the  local  geometry  from  which  each  vortex 
was  generated.  The  different  shedding  characteristics  might  lead  to  the  mutual  dissipation  of 
coherent  vortical  structures,  which  eventually  might  result  in  higher  flame  stability. 

Along  the  x-coordinate,  two  vortical  structures,  ux\  and  uix 2,  also  were  formed:  by  the  back 
steps  with  the  side  surface  of  the  S  cavity  and  the  top  surface  of  the  no-cavity  region  (wTl), 
and  the  side  surface  of  the  S  cavity  and  the  bottom  surface  of  the  D  cavity  (lox 2),  respectively. 
We  believed  that  these  vortices,  along  with  the  uz 2,  were  essential  to  mixing  between  fuel  that 
flowed  through  the  S  cavity  region  and  air,  which  flowed  along  the  no-cavity  region  (u+1)  or 
resided  inside  of  the  D  cavity  ( ux2  and  ujz2). 

Finally,  a  vortical  structure  along  the  y-coordinate,  ujy  1,  was  also  present,  as  previously 
confirmed  in  Figs.  10  and  14.  This  vortex  was  produced  by  an  inter-cavity  flow  that  started 
from  one  of  the  S  cavities,  passed  through  the  inside  of  the  D  cavity,  and  finally  exited  to  the 
left  corner  of  the  adjacent  cavity  on  its  right.  As  a  result,  the  flow  produced  a  swirling  motion 
surrounding  the  no-cavity  region.  The  swirling  motion  was  clockwise  when  seen  from  the  top, 
and,  therefore,  the  resulting  vortex  pointed  in  the  — y  direction.  We  believed  that  its  companion 
vortex  (pointing  in  the  +y  direction)  should  be  present  if  ideal  symmetry  could  be  achieved. 

For  the  T+S  geometry,  a  large  vortex  due  to  the  slanted  side  surface  of  the  T  geometry 
(ujz  1),  a  small  vortex  due  to  the  back  step  formed  by  the  side  surface  of  the  bluff  body  and 
the  bottom  surface  of  the  S  cavity  (<n23),  and  a  mid-size  vortex  due  to  the  back  step  formed  by 
the  back  surface  of  the  S  cavity  and  the  slanted  T  surface  (c+2)  were  present.  An  additional 
vortex  (ux  1)  was  present  in  this  geometry  along  the  x  coordinate  due  to  the  back  step  formed 
by  the  side  surfaces  of  the  S  cavity  and  the  slanted  T  surface.  This  vortex  induced  inter-cavity 
flow  along  the  slanted  T  surface  observed  in  Figs.  14  and  15.  Again,  these  multiple  vortices  had 
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(a) 


Figure  18:  Overall  flow  pattern  of  (a)  the  Reference,  (b) 


the  D+S  and  (c)  the  T+S  geometries. 
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different  shedding  frequencies  and  locations. 


2.3  Blowout  Study 

Here,  we  present  variations  of  the  flame  configuration  to  study  blowout  effects.  In  this  section, 
our  discussion  will  be  limited  to  the  six  base  geometries:  R,  D,  T,  R+S,  D+S,  and  T+S  shown 
in  Fig.  4.  The  results  from  these  selected  geometries  would  show  the  significance  of  the  S  cavity. 

2.3.1  Variation  of  Flame  Configuration 

Three  flame  configurations  were  considered:  i)  a  fully  premixed  freestream  stabilized  by  a  bluff 
body  (fully  premixed  mode);  ii)  an  air  freestream  in  which  fuel  was  discharged  from  the  bluff 
body  (partially  premixed  mode);  and  iii)  a  combination  of  these  two,  namely  a  fully  premixed 
freestream  with  additional  fuel  discharged  from  the  bluff  body.  The  premixed  burner  that  was 
used  to  provide  the  vitiated  air  flow  in  the  liftoff  study  also  was  used  here  to  supply  a  fully 
premixed  (for  fully  premixed  and  hybrid  modes)  or  pure  air  (for  partially  premixed  mode) 
stream  at  room  temperature. 

2.3.2  Gas  Chromatography  for  Measuring  Methane  Concentration 

Gas  chromatographic  measurement  was  performed  to  quantify  the  local  methane  concentration 
in  the  vicinity  of  the  bluff  body  base  using  a  Varian  3400  Gas  Chromatograph  equipped  with 
Porapak  Q  (for  nitrogen  and  oxygen)  and  Molecular  Sieve  5A  (for  methane)  as  columns  and  a 
thermal  conductivity  detector  (TCD).  All  sampling  measurements  were  conducted  in  the  absence 
of  the  flame.  For  the  sampling  the  spatial  scanning  was  carried  out  by  three  translation  stages 
(along  the  x,  y,  and  z  coordinates)  using  a  ~1  mm  diameter  (ID),  200  mm  long  stainless  steel 
probe. 

2.3.3  Blowout  Limit  Measurements 

Figure  19  shows  the  results  obtained  from  the  ensemble  averaged  flame  blowout  limit.  These 
results  were  quantified  in  terms  of  the  local  flow  speed  measured  at  2  cm  distance  away  from 
the  methane  nozzle  for  all  six  cases  discussed  in  Section  2.3.1  operating  in  the  hybrid,  fully 
premixed,  and  partially  premixed  conditions.  To  vary  the  local  flow  speed,  additional  air  from 
the  turbulent  grid  was  injected,  and  its  flowrate  was  varied  (from  40  SLPM  to  200  SLPM)  until 
the  entire  bluff  body  stabilized  flame,  which  was  formed  from  the  fully  premixed  flow  from 
the  gas  blower  and  pure  methane  jet  injected  from  the  methane  nozzle,  was  blown  out.  The 
equivalence  ratio  and  the  speed  of  the  fully  premixed  stream  from  the  gas  blower  were  kept 
identical  at  0.72  and  ~0.83  m/s,  respectively.  The  amount  of  methane  injection  from  the  nozzle 
also  was  kept  fixed  at  ~3  SLPM  such  that  the  (nominal)  square  root  of  the  momentum  ratio 
(r)  calculated  from  the  main  fully  premixed  stream  and  the  injected  methane  was  ~1.2  (the  r 
parameter  ranges  from  0.2  to  0.9  when  it  was  calculated  based  on  the  nominal  local  speed  formed 
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by  the  main  fully  premixed  flow  and  additional  air  injection).  For  the  blowout  experiments,  the 
relative  placement  of  the  methane  nozzle  and  the  local  cavity  was  in-phase. 

The  critical  nominal  flow  speeds  in  the  blowout  limit  (black  solid  bar  in  Fig.  19a)  were  as 
low  as  ~3.7  m/s  for  R,  ~3.7  m/s  for  D,  and  ~4.7  m/s  for  T  geometies,  which  did  not  include 
any  cavity.  The  critical  speed  was  improved  significantly  by  adding  local  cavities  to  the  original 
geometries;  i.e.,  by  48%  for  the  R+S  geometry,  by  42%  for  the  D+S  geometry,  and  by  10% 
for  the  T+S  geometry.  The  difference  between  the  T  and  T+S  geometries  was  small,  which 
probably  was  due  to  the  fact  that  i)  the  T  geometry  already  showed  a  good  capability  to  extend 
the  blowout  limit  in  comparison  with  those  of  R  and  D  and  ii)  the  size  of  the  S  cavity  on  the 
T  geometry  was  small  (due  to  slanted  side  surface  of  T),  and,  therefore,  the  effect  of  the  S 
cavity  was  relatively  weaker.  The  identical  experiments  without  the  injection  of  pure  methane 
(fully  premixed  mode),  denoted  by  the  shaded  bars  in  Fig.  19a,  showed  no  apparent  stability 
enhancement  with  addition  of  the  S  cavity.  This  observation,  i.e.,  minimal  positive  effect  on 
blowout  limit  extension  with  the  geometrical  change  of  bluff  body,  was  consistent  with  previously 
published  results  (Stwalley  &  Lefebvre  (1988)  and  Khosla  et  al.  (2007)).  For  comparison,  we 
also  carried  out  an  experiment  that  is  similar,  but  with  methane  injection  into  pure  air  crossflow 
(partially  premixed  mode).  As  shown  in  the  dotted  bar,  no  significant  extension  of  the  blowout 
limit  was  observed  among  the  various  geometries. 


R  D  T  R+S  D+S  T+S  R  D  T  R+S  D+S  T+S 


Figure  19:  Blowout  limits  in  terms  of  (a)  local  flow  speed  (measured  at  2  cm  horizontal  distance 
from  the  methane  nozzle)  and  (b)  equivalence  ratio  of  fully  premixed  stream  for  various  bluff 
body  base  designs.  For  (a),  the  equivalence  ratio  of  the  fully  premixed  stream  was  kept  fixed  as 
0.72  (for  the  fully  premixed  and  hybrid  modes)  whereas  in  (b)  nominal  local  flow  speeds  of  were 
kept  identical  to  1.5  m/s  (fully  premixed  mode)  and  2.8  m/s  (hybrid  mode). 


In  Fig.  19b,  the  blowout  limits  in  terms  of  equivalence  ratio  (of  the  fully  premixed  main 
stream)  for  all  six  geometries  are  plotted.  The  injected  air  rates  (from  the  turbulent  grid)  were 
kept  fixed  at  ~50  SLPM  (in  the  fully  premixed  case)  and  ~100  SLPM  (in  the  hybrid  case), 
respectively.  The  blowout  limit  was  measured  by  gradually  decreasing  the  injected  methane 
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rate  in  the  fully  premixed  stream  until  no  bluff  body-stabilized  flame  was  present  while  the  air 
flowrate  of  the  stream  was  kept  identical  at  ~1720  SLPM.  The  resulting  nominal  flow  speeds 
of  the  stream  were  1.5  m/s  and  2.8  m/s  for  the  fully  premixed  mode  and  the  hybrid  mode, 
respectively. 

In  the  hybrid  mode  (black  solid  bar),  the  presence  of  the  local  cavity  in  the  R+S,  D+S,  and 
T+S  geometries  extended  the  blowout  limit  from  0.62  (R  and  D)  and  0.61  (T)  to  0.55  (R+S), 
0.56  (D+S)  and  0.58  (T+S).  This  trend  was  consistent  with  what  was  observed  in  Fig.  19a. 
Again,  for  the  fully  premixed  mode  (shaded  bar),  no  significant  dependence  of  the  blowout 
limits  on  the  geometries  tested  was  observed,  which  also  confirmed  the  previous  observation 
in  Fig.  19a.  Thus  Figs.  19a  and  19b  suggested  that  there  were  no  clear  benefits  to  blowout 
limits  in  both  the  fully  premixed  and  partially  premixed  modes  of  operation,  but  there  were 
clear  benefits  in  the  hybrid  mode  of  operation.  The  fully  premixed  results  were  consistent  with 
observations  of  Stwalley  &;  Lefebvre  (1988)  and  might  suggest  that  the  presence  of  large-scale 
oscillations  was  not  an  indication  of  blowout  and  that  this  phenomenon  was  more  likely  related 
to  the  observations  of  Shin  &  Ferziger  (1991). 

2.3.4  Measurements  of  Methane  Mole  Fraction  Field 

Figures  20a  and  20b  show  the  variation  of  methane  mole  fraction  along  the  y  and  z  coordinates 
in  the  non-burning  hybrid  mode  case,  respectively.  For  all  of  these  mole  fraction  measurements, 
the  equivalence  ratio  and  the  speed  of  the  fully  premixed  stream  were  kept  fixed  as  0.72  and 
0.83  nr/s,  respectively.  The  flowrate  of  the  additional  air  injection  also  was  kept  constant  at 
100  SLPM.  The  resulting  speed,  including  the  air  injection  (measured  at  x  «20  mm),  was  ~3 
nr/s.  The  measurements  in  Fig.  20a  were  taken  at  x  =  0  mm,  the  center  of  bluff  body  in  the  xy 
plane,  and  z  =  0  mm,  the  center  of  local  cavity  in  the  yz  plane.  The  measurements  in  Fig.  20b 
are  taken  at  x  =  8  mm  and  y  =  6.5  mm,  from  the  top  surface  of  the  bluff  body. 

Along  the  y  coordinate  (Fig.  20a),  the  methane  mole  fractions  of  the  non-local  cavity  geome¬ 
tries,  especially  R  and  D,  were  generally  lower  by  ~50  %  compared  to  those  of  the  local  cavity 
geometries  (R+S,  D+S  and  T+S)  at  least  up  to  y  =  ID,  where  D  was  the  width  of  the  bluff  body 
(1  inch).  We  believed  that  the  high  fuel  mole  fractions,  especially  in  the  vicinity  of  the  base, 
played  an  important  role  in  increasing  the  flame  stability  for  the  local  cavity  geometries.  We 
expected  that  the  most  active  interactions  between  the  base  geometry  and  the  flow  would  occur 
in  the  vicinity  of  the  base.  The  results  showed  that  the  presence  of  the  local  cavity  prevented 
the  pure  methane  stream  from  being  diluted  excessively  for  which  the  partially  premixed  stream 
effectively  could  pilot  the  fully  premixed  stream.  The  T  geometry  showed  a  mole  fraction  field 
as  high  as  those  of  local  cavity  geometries.  This  high  mole  fraction  was  clear  in  the  PIV  results 
(see  Fig.  21  for  the  detailed  flow  field),  which  indicated  the  presence  of  a  skewed  flow  along  the 
slanted  side  surface  of  the  T  geometry.  We  believed  that  the  skewed  flow  effectively  could  deliver 
the  pure  methane  stream,  especially  at  this  low  value  of  r  (r  <  1),  to  the  region  investigated 
here  such  that  the  overall  fuel  mixture  fraction  in  that  region  remained  high. 
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Figure  20:  Methane  mole  fractions  along  the  (a)  y  and  (b)  z  directions  for  the  six  modified 
geometries.  The  equivalence  ratio  of  the  fully  premixed  stream  was  kept  fixed  as  0.72.  The 
flowrate  of  the  additional  air  injection  also  was  kept  at  100  SLPM  such  that  the  resulting 
nominal  flow  speed  measured  at  2  cm  distance  from  methane  nozzle  was  ~3  m/s.  Measurement 
locations  were  x  =  0  and  z  =  0  for  (a)  and  x  =  8  mm  and  y  =  6.5  mm  for  (b). 


Figure  20b  shows  results  of  the  same  experiment,  as  described  in  the  preceding  paragraph, 
but  along  the  z  coordinate.  From  the  figure,  a  similar  trend  existed,  i.e. ,  the  local  cavity 
geometries  (and  the  T  geometry)  contained  approximately  two  times  higher  local  fuel  mole 
fraction  compared  with  the  non-local  cavity  geometries  (especially  R  and  D).  The  mole  fractions 
of  R+S  (black  dotted  line)  and  D+S  (red  dotted  line)  geometries  decreased  with  increasing  z, 
while  the  other  four  geometries  showed  a  relatively  invariant  tendency  of  mole  fraction  along 
the  z  direction.  The  decay  was  more  apparent  when  z  was  greater  than  ~6  mm.  This  decay  was 
because  the  local  cavity  section  ended,  and  the  wall  section  started  from  z  =  6.4  mm  in  those 
two  geometries,  which  confirmed  that  the  local  cavity  acted  as  an  effective  conduit  of  less  diluted 
fuels.  To  study  the  capability  of  local  cavity  placement  relative  to  the  methane  jet  to  localize  the 
fuel  stream  further,  we  briefly  investigated  the  blowout  limit  dependence  (under  experimental 
conditions  of  Fig.  20)  on  the  location  of  the  cavity  relative  to  the  methane  jet  nozzles.  The 
result  showed  that  the  blowout  limit  for  out-of-phase  placement  (i.e.  the  case  that  the  center 
of  the  wall  section  was  aligned  with  the  methane  nozzle)  was  ~20  %  (in  terms  of  the  local  flow 
speed)  or  ~7%  (in  terms  of  the  equivalence  ratio  of  fully  premixed  stream)  worse  than  those 
of  in  phase  placement  (i.e.,  the  case  where  the  center  of  the  cavity  section  was  aligned  with  the 
methane  nozzle)  for  both  the  R+S  and  D+S  geometries. 

It  was  important  to  compare  the  mole  fraction  and  velocity  fields  simultaneously.  The 
total  amount  of  fuel  injection  was  identical  for  all  the  studied  geometries,  and,  therefore,  the 
lower  methane  mole  fraction  for  non-local  cavity  geometries  implied  that  there  were  regions  of 
higher  mole  fraction  (see  Fig.  20).  Hence,  we  believed  that  the  distribution  of  the  fuel  mole 
fraction  in  the  flow  held  was  likely  more  important  than  the  absolute  value  of  the  fuel  mole 
fraction  (shown  in  Fig.  20)  itself.  To  investigate  the  mole  fraction  distribution  resulting  from 
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Figure  21:  (a)  Instantaneous  and  (b)  ensemble  averaged  PIV  images  of  the  R,  D,  T,  R+S,  D+S, 
and  T+S  geometries.  The  nominal  flow  speed  was  ~3  m/s  and  the  equivalence  ratio  of  the 
fully  premixed  stream  was  kept  fixed  as  0.72  for  all  images.  The  red  lines  represent  the  pseudo 
streamlines. 
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different  flow  fields,  we  carried  out  gas  chromatographic  sampling  in  parallel  to  independent  PIV 
measurements.  The  gas  chromatographic  sampling  results  are  shown  in  Fig.  22.  In  this  figure, 
the  variations  of  methane  mole  fraction  along  the  x  coordinate  (at  y  =  6.5  mm  and  z  =  0  mm) 
were  plotted  together  with  velocity  components,  Vx  and  Vy,  at  the  equivalent  locations  of  the 
gas  sampling  for  all  the  six  geometries.  Similar  to  the  case  in  Fig.  20,  both  the  sampling  and 
PIV  measurements  were  for  the  non-reacting  case. 

As  shown  in  Fig.  22,  for  the  R  (Fig.  22a)  and  D  geometries  (Fig.  22b),  fuel  mole  fractions 
were  relatively  low  (~1.8  -  3.3  %)  at  x  <~  10  mm,  where  the  Vx  and  Vy  velocities  were  also  very 
low  (<~0.1  m/s),  while  the  higher  velocity  (Vy  ~  3  m/s)  region  showed  higher  fuel  mole  fraction 
(~5.5  %).  This  type  of  fuel  distribution,  i.e.  high  fuel  concentration  in  a  high  velocity  region  and 
low  fuel  concentration  in  a  low  velocity  region,  was  not  the  most  efficient  way  to  improve  flame 
stability.  To  the  contrary,  for  the  local  cavity  geometries  in  Figs.  22d  -  22f  (R+S,  D+S,  and 
T+S,  respectively),  the  high  fuel  mole  fraction  («  5  -  7%)  region  was  maintained  in  relatively 
low  velocity  regions  (Vy  ~  1  —  2  m/s,  Vx  ~  0  -  1  m/s),  where  x  <«10  mm,  whereas  the  mole 
fraction  decreased  in  the  high  velocity  region  (x  >~15  mm).  For  example,  the  region  of  peak 
mole  fraction  in  the  R+S  geometry  (x  ~4  mm)  contained  regions  with  velocities  Vy  ~  2  m/s 
(69%  of  peak  Vy)  and  Vx  ~  0.8  m/s,  whereas  in  the  R  geometry  (x  ~16  mm)  the  velocities 
were  Vy  ~  3  m/s  (99%  of  peak  Vy)  and  Vx  ss  0.1  m/s. 

An  interesting  observation  could  be  made  by  looking  at  the  gradients  of  fuel  mole  fraction. 
For  the  R  and  D  geometries,  the  mole  fractions  increased  by  more  than  100  %  (from  ~2.5  %  to 
~6  %)  within  ~4-mm  distance,  but  for  the  local  cavity  geometries,  the  variations  of  the  mole 
fraction  were  much  smaller,  less  than  25  %  in  the  highest  case.  We  believed  that  specific  fuel 
distributions,  as  well  as  smaller  mole  fraction  gradient  in  the  presence  of  the  local  cavities, 
effectively  improved  blowout  limits  in  the  hybrid  mode  flow. 

As  a  final  comment  on  the  results  presented  in  Fig.  22,  we  compared  the  mole  fraction  and 
velocity  fields  of  T  (Fig.  22c)  and  T+S  (Fig.  22f)  geometries.  The  behavior  of  the  two  fields  was 
generally  similar  for  both  geometries.  Both  cases  showed  a  gradual  decay  of  the  mole  fraction 
field  from  «6  %  to  4.5  %,  a  gentle  increase  of  Vy  up  to  x  ~12.5  mm,  and  an  asymptotic  decrease 
to  the  free  stream  velocity  with  increasing  x.  As  mentioned  earlier,  this  similarity  was  caused 
by  the  fact  that  the  size  of  the  cavity  in  the  T+S  geometry  was  relatively  small;  however,  one 
could  still  observe  some  evidence  of  the  presence  of  the  local  cavity  from  the  mole  fraction  data 
at  x  ~12  mm  in  the  T+S  geometry  (orange  circle).  The  mole  fraction  of  the  T+S  geometry 
was  10%  higher  than  that  of  T,  while  the  neighboring  mole  fractions  (e.g.,  at  x  ~10  mm  and 
14  mm)  were  almost  identical  for  both  geometries.  We  believed  that  the  presented  data  further 
confirmed  the  capability  of  a  local  cavity  to  localize  the  fuel  stream. 


34 


X  (mm) 


* 


It 


7 

£  6 
c 

|  5 
u 
2 

Q)  4 


4  i 


°  0 
0* 

> 

-1 


X  (mm) 


X  (mm) 


X  (mm) 


Figure  22:  Methane  mole  fractions  and  velocity  components  (Vx  and  Vy)  along  the  x  direction 
for  six  geometries:  (a)  R,  (b)  D,  (c)  T,  (d)  R+S,  (e)  D+S  and  (f)  T+S.  Experimental  conditions 
were  identical  to  those  in  Fig.  20.  Measurement  locations  were  at  y  =  6.5  and  z  =  0. 
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2.3.5  Flow  Fields  Near  the  Base  in  a  Room  Temperature  Freestream  (2D  PIV 
Measurement ) 

Figure  21  shows  instantaneous  (Fig.  21a)  and  ensemble  averaged  (Fig.  21b)  PIV  images  for  the 
six  non-cavity  and  cavity  geometries,  i.e.,  R,  D  and  T,  R+S,  D+S,  and  T+S.  Also  shown  in 
red  are  pseudo-streamlines.  For  all  the  cases  shown  in  Figure  21,  the  flame  was  present  and  the 
configuration  corresponded  to  the  hybrid  mode.  The  flow  rate  of  additional  air  injection  and  the 
resulting  nominal  flow  speed  (measured  at  2  cm  distance  away  from  methane  nozzle)  were  kept 
identical  as  100  SPLM  and  ~3  m/s,  respectively.  The  equivalence  ratio  of  the  fully  premixed 
stream  also  was  kept  fixed  at  0.72  for  all  images. 

As  shown  in  the  instantaneous  images  (especially  for  R,  D,  R+S,  and  D+S  geometries), 
recirculating  flow  structures,  which  were  believed  to  play  a  key  role  in  the  bluff  body  flame 
stabilization,  could  be  found  for  all  four  geometries.  For  example,  one  could  see  clearly  strong 
recirculating  pairs  with  their  centers  located  at  y  ~10  mm  and  y  ~6  mm  in  the  R  geometry  and 
at  y  ~10  mm  and  y  ~20  mm  in  the  D  geometry.  For  the  R+S  and  D+S  geometries  the  structure 
resided  at  the  height  of  y  ~5  mm  (R+S)  and  y  ~-2  mm  (D+S).  The  recirculating  pattern  was 
more  clear,  and  its  strength  was  more  intense  in  the  non-cavity  geometries  than  in  the  cavity 
cases.  The  recirculation  in  the  cavity  cases  was  less  visible  because  either  intense  recirculation 
zones  of  the  local  cavity  geometries  resided  in  a  region  where  the  laser  light  could  not  illuminate, 
e.g.,  inside  of  the  cavity,  or  the  geometries  tended  to  generate  small,  multiple  vortices  rather 
than  one  large  vortex  pair  (see  Fig.  23  for  a  more  detailed  discussion  on  the  characteristics  of 
vortices  generated  by  each  geometry). 

The  vortices  generated  from  the  local  cavity  geometries  were  more  incoherent,  i.e.,  the  shed¬ 
ding  frequency,  size,  and  locations  were  more  diverse,  as  described  in  the  liftoff  study.  Evidence 
of  the  incoherent  vortices  also  could  be  seen  in  the  present  study  as  shown  in  the  averaged  PIV 
images  of  Fig.  21b.  The  recirculation  pattern  for  non-cavity  geometries  was  still  clear  from  this 
figure  (especially  for  R  and  D  geometries);  however,  no  apparent  vortex  was  observed  in  the 
cavity  geometries  (R+S  and  D+S).  This  fact  implied  that  the  vortices  observed  in  the  instanta¬ 
neous  images  of  the  cavity  geometries  were  more  random  (incoherent)  and  therefore  tended  to 
disappear  in  the  ensemble- averaged  images. 

Unlike  the  four  geometries  D,  R,  D+S,  and  R+S,  the  T  and  T+S  geometries  had  more  direct 
(i.e.,  non-recirculating)  flow  patterns  along  the  slanted  side  surfaces  in  the  region  of  interest  for 
which  laser  illumination  was  present.  The  flow  pattern  was  inconsistent  with  our  previous 
observations  (liftoff  study),  which  showed  two  large,  nearly  symmetric  recirculation  zones  along 
the  side  surfaces  for  the  ~5.7  m/s  vitiated  flow  in  the  absence  of  additional  air  injection.  Slight 
misalignment  between  the  direction  of  fully  premixed  flow  and  the  air  injection  in  the  current 
configuration  induced  a  small  skew  of  the  local  flow  to  the  right  in  the  region  of  interest.  The 
skewed  flow,  as  clearly  seen  in  the  averaged  field,  seemed  to  eliminate  the  recirculation  zone  on 
the  left  face  in  the  region  of  interest  and  created  bent  flow  along  the  vortices  at  y  ~0  (clockwise 
direction).  There  was  evidence  of  a  recirculation  zone  on  the  right  face  of  the  T  and  T+S  bases 


36 


10 


Figure  23:  Swirling  strength  of  (a)  instantaneous  and  (b)  ensemble  averaged  velocity  fields  of 
the  R,  D,  T,  R+S,  D+S,  and  T+S  geometries  obtained  under  experimental  conditions  identical 
to  Fig.  21.  The  white  lines  are  pseudo  streamlines. 


(e.g.,  near  x  =  y  =  -5)  from  the  relatively  upright  averaged  flow  field  at  x  ~-5  -  -10,  y  =  0 
in  Fig.  21b,  and  the  corresponding  instantaneous  images  in  Fig.  21a.  This  small  asymmetry 
might  be  amplified  in  geometries  such  as  T  and  T+S  (which  showed  a  higher  sensitivity  to 
the  flow  symmetry  because  of  their  sharp  vortices)  and  thus  altered  the  flow  pattern  of  those 
geometries,  resulting  in  either  two  or  just  a  single  recirculation  bubble.  This  aspect  did  not 
seem  to  change  our  key  results  in  the  current  study  (i.e.,  improvement  of  the  blowout  limit  in 
the  cavity  geometries  due  to  geometry’s  ability  to  distribute  the  fuel  in  a  less  diluted  manner) 
in  that  our  observations  still  held  for  geometries  such  as  the  R  and  R+S  and  the  D  and  D+S, 
which  were  more  robust  to  the  flow  asymmetry. 

2.3.6  Swirling  Strengths  and  Strain  Rates  Near  the  Base  in  a  Room  Temperature 
Freestream 

Figure  23a  shows  instantaneous  images  of  swirling  strength  field  (and  pseudo  streamlines  rep¬ 
resented  by  white  lines)  for  the  six  geometries  R,  D,  T,  R+S,  D+S,  and  T+S  at  experimental 
conditions  identical  to  those  of  Fig.  21.  It  was  known  that  the  swirling  strength,  a  measure  of 
how  fast  the  fluid  is  rotating  locally,  was  a  good  marker  of  vortical  structures  (see  Christensen 
&  Wu,  2005).  There  were  more  vortices  in  the  bluff  body  downstream  of  the  cavity  geometry 
(and  T  geometry)  than  in  non-cavity  geometries  (especially  in  the  R,  R+S,  D,  and  D+S).  We 
believed  that  more  complex  back-steps  in  the  geometry  helped  generate  multiple  and  incoherent 
vortices. 

In  Fig.  23b  the  swirling  strength  of  the  averaged  velocity  field  (shown  in  Fig.  21b)  is  il¬ 
lustrated,  where  the  swirling  strength  was  defined  as  the  imaginary  portion  of  the  complex- 
conjugate  eigenvalues  of  the  local  velocity  gradient  tensor  (see  Natrajan  &  Christensen,  2006). 
Unlike  the  instantaneous  field  in  Fig.  23a,  no  noticeable  discrepancy  existed  between  the  non¬ 
cavity  and  cavity  geometries,  i.e.,  there  existed  a  similar  number  of  diffuse  vortices  in  all  six 
cases.  This  observation  (multiple,  highly  concentrated  vortices  in  an  instantaneous  field  and 
a  smaller  number  of  diffuse  vortices  in  an  averaged  field  shown  in  the  cavity  geometries)  con¬ 
firmed  that  the  local  cavity  geometry  could  generate  incoherent  vortices  unlike  the  non-cavity 
geometries  (especially  in  the  R  and  D  geometries). 

Finally,  the  strain  rates  of  the  averaged  velocity  field  in  the  six  representative  geometries  are 
illustrated  in  Fig.  24.  The  experimental  conditions  were  identical  to  those  of  Fig.  23.  As  shown 
in  this  figure,  the  regions  of  high  strain  rates  were  generally  large  in  the  non-cavity  geometry 
cases  (e.g.,  the  strain  rate  on  the  top  surface  of  the  R  and  D  geometries),  whereas  that  region 
was  small  in  cavity  geometry  (R+S  and  D+S).  We  did  not  believe  that  this  level  of  strain  rate 
(<~100/ s)  could  cause  local  extinction,  but  we  expected  that  the  high  strain  rate  in  the  vicinity 
of  the  bluff  body  seen  in  the  non-cavity  geometries  did  not  increase  the  flame  stability.  Similar 
to  previous  observations,  no  noticeable  difference  was  found  between  the  T  and  T+S  geometries. 
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Figure  24:  Strain  rate  obtained  from  ensemble  averaged  velocity  field.  Experimental  conditions 
were  identical  to  that  of  Fig.  21.  The  white  lines  represent  the  pseudo  streamlines. 
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3  Numerical  Modeling  Using  Large-Eddy  Simulation 


In  the  previous  section,  several  new  flame  holder  designs  were  presented  with  improved  stability 
characteristics.  The  experimental  characterization  of  the  developed  geometries  was  done  in  a 
low-speed  vitiated  flow  configuration.  In  the  present  section,  two  different  approaches  were 
presented  to  assist  in  the  development  of  new  flame-holding  technologies.  The  first  was  a  new 
formulation  of  Damkbhler  (Da)  number  similarity.  This  technique  allowed  for  a  proper  definition 
of  Da,  such  that  similarity  conditions  could  be  defined  for  the  assessment  of  flame  stabilization. 
With  this  method,  low-speed  experiments  could  be  used  to  assess  afterburner  static  stability 
under  typical  operating  conditions.  The  second  was  an  assessment  of  the  combustion  modeling 
techniques  to  predict  the  static  stability  of  the  flame-holder  configurations  presented  in  Section 
2. 

3.1  Damkohler  Number  Similarity  for  Static  Flame  Stability  in  Augmentor 
Flows 

3.1.1  Introduction 

Flame  holders  provided  several  features  important  for  flame  stability.  As  the  presence  of  the 
flame  holder  modified  the  velocity  field,  a  recirculation  zone  (RZ)  was  generated  after  the  flame- 
holder  base.  This  recirculation  zone  provided  sufficient  residence  time  and  chemical  heat  release 
needed  to  ignite  the  fresh  gas  mixture.  In  addition,  the  shear  layer  generated  around  the  RZ 
provided  a  turbulent  mechanism  to  mix  fresh  fuel  and  air  with  the  combustion  products  from 
the  recirculation  zone  (Zukowski  &;  Marble  (1955);  Longwell  et  al.  (1953)).  Finally,  as  the  flame 
was  stabilized  past  the  flame  holder  base,  combustion  products  at  the  end  of  the  recirculation 
zone  were  entrained  into  the  RZ  to  provide  a  self-ignition  mechanism. 

The  inflow  mixture  to  the  augmentor  usually  was  at  high-speed  and  high-temperature.  The 
typical  inflow  temperatures  and  velocities  for  augmentors  ranged  between  700  and  1000  K  and 
100  and  200  m/s,  respectively  (King  &  Nakanishi  (1957)).  Under  these  conditions  the  experimen¬ 
tal  assessment  of  such  configurations  was  complex  and  expensive.  The  current  work  investigated 
the  possibility  of  utilizing  mathematical  similarity  to  test  such  configurations.  Here,  we  assessed 
the  Damkohler  number  similarity,  which  implied  that,  in  the  high  Reynolds  number  limit,  the 
flow  and  the  flame  characteristics  in  a  high-speed/fast-chemistry  case  should  be  similar  to  a 
low-speed/slow-chemistry  case,  as  long  as  the  Damkohler  number  and  all  the  momentum  ratios 
were  the  same. 

Several  papers  have  quantified  flame  blow-off,  and  the  Damkohler  number  often  was  used  in 
the  parameterization  (Mellor  (1980);  Radhakrishnan  et  al.  (1981)).  A  comprehensive  review  of 
the  subject  and  a  collection  of  data  for  bluff-body-stabilized  flames  in  augmentor-like  systems 
recently  was  provided  by  Shanbhogue  et  al.  (2009a).  An  early  attempt  to  formulate  a  general 
theory  of  premixed  flame  stabilization  behind  a  bluff  body  was  made  by  Longwell  et  al.  (1953), 
who  derived  flame  stability  curves  for  different  activation  energies  assuming  a  second-order 
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homogeneous  reaction  at  steady-state  conditions.  King  &  Nakanishi  (1957)  investigated  the 
effect  of  the  flame  holder  geometry  on  combustion  efficiency.  They  reported  a  strong  correlation 
between  fuel  to  air  ratio  at  extinction  and  the  inflow  parameters  of  pressure,  temperature, 
and  velocity.  More  recently,  two  experimental  and  numerical  studies  (Kiel  et  al.  (2007,  2006)) 
discussed  flame  extinction  behind  a  set  of  bluff-body-stabilized  flames.  They  speculated  that 
three  different  effects  controlled  wake-stabilized  flames.  These  effects  were  viscous  forces,  inertial 
forces,  and  the  baroclinic  torque  generated  by  the  density  variation  due  to  combustion.  The 
competition  among  these  effects  changed  the  flame  blow-out  and  extinction  characteristics.  For 
high  Reynolds  numbers,  in  the  square,  V-gutter,  and  circular  cylinder  flame  holders,  the  latter 
two  effects  were  found  to  be  most  important.  At  low  Reynolds  numbers,  the  viscous  forces  were 
found  to  be  more  dominant. 

The  near-blow-off  dynamics  of  a  bluff-body-stabilized  flame  were  studied  by  Lieuwen  et  al. 
(2007)  and  Suraj  &  Lieuwen  (2007).  The  unsteadiness  was  found  to  be  a  major  factor  that 
promotes  blow-off.  Lieuwen  and  Suraj  concluded  that  the  flame  blowout  could  not  be  explained 
based  on  purely  steady  state  characteristics  such  as  the  residence  and  the  chemical  time  scales. 
They  argued  that  two  factors  mainly  drove  flame  extinction,  the  presence  of  vortex  shedding  in 
the  wake  of  the  bluff-body  and  the  locally  induced  strain.  The  combined  effects  of  these  two 
unsteady  factors  controlled  the  path  of  the  reactants  and  the  hot  gases,  and,  under  certain  con¬ 
ditions,  if  the  recirculation  was  strong  enough,  the  cold  gases  could  sweep  the  stabilization  point, 
leading  to  flame  blow-off.  In  that  sense,  the  onset  of  the  large-scale  deterministic  unsteadiness 
could  be  seen  as  a  precursor  for  blow-off;  therefore,  in  this  report,  although  the  main  interest 
was  in  actual  blow-off,  we  only  considered  the  near-blow-off  regime,  which  was  characterized  by 
vortex  shedding.  We  then  used  the  characteristics  of  that  large-scale  unsteadiness  as  a  surrogate 
in  the  discussion  of  static  flame  stability  without  determining  the  actual  stability  limits. 

Several  computational  studies  of  bluff-body  flame  stabilization  have  been  reported.  Porumbel 
&  Menon  (2006),  for  instance,  performed  an  LES  study  of  a  bluff-body-stabilized  premixed 
flame.  They  concluded  that  the  baroclinic  torque  generated  by  combustion  weakened  the  Kelvin- 
Helmholtz  instability  and  suppressed  the  von  Karrnan  vortex  street  observed  in  the  non-reactive 
flow  simulation.  The  same  conclusion  was  supported  by  Mehta  &  Soteriou  (2003)  and  Ghoniem 
&  Krishnan  (1988).  Khosla  et  al.  (2007)  studied  a  bluff-body  configuration  similar  to  the  current 
work.  They  highlighted  the  effect  of  suppressing  the  von  Karrnan  shedding  on  the  flame  blowout 
characteristics.  By  using  a  tabbed  bluff-body,  the  coherence  of  the  vortex  shedding  between  the 
upper  and  the  lower  surface  of  the  bluff-body  was  disrupted,  such  that  the  asymmetrical  vortex 
shedding  mode  could  be  eliminated.  Knaus  et  al.  (2008)  and  Roach  et  al.  (2008)  used  different 
numerical  approaches  in  FLUENT  ( www.ansys.com )  to  predict  the  Damkohler  number  in  a  V- 
gutter-stabilized  flame.  They  defined  the  local  chemical  time  scale  based  on  the  extinction  strain 
rate  using  a  single-step  global  chemistry.  This  extinction  strain  rate  was  found  to  increase  with 
the  global  equivalence  ratio. 

All  these  experimental  and  computational  studies  have  contributed  to  the  understanding  of 
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bluff-body  supported  flame  stabilization;  however,  there  still  is  a  need  for  fundamental  studies 
and  flame  holder  testing.  Experimental  investigations  can  be  simplified  greatly  by  considering 
Damkohler  similarity.  The  basic  idea  is  that  all  flow  and  combustion  characteristics,  including 
stability,  should  be  similar  for  different  operating  conditions,  as  long  as  all  non-dimensional 
groups  describing  the  problem  are  the  same.  One  of  the  more  important  non-dimensional  groups 
important  for  combustion  problems  is  the  Damkohler  number,  which  is  defined  as  the  ratio  of 
characteristic  flow  and  chemistry  time  scales;  however,  changing  the  operating  conditions  while 
keeping  the  Damkohler  number  constant  is  not  trivial.  In  this  report,  a  method  for  ensuring 
Damkohler  number  similarity  was  provided  and  validated. 

Numerical  simulations  for  three  different  conditions  were  performed  using  large-eddy  sim¬ 
ulations,  where  combustion  was  modeled  using  the  flamelet  progress  variable  (FPV)  approach 
(Pierce  &;  Moin,  2004).  The  configuration  was  for  a  bluff-body-stabilized  flame  in  a  vitiated  flow. 
Methane  fuel  was  injected  in  the  cross-stream  direction  from  the  bluff-body  surface.  Simulation 
results  for  a  reference  case  were  compared  with  experimental  data,  which  also  were  presented  in 
this  report.  The  inflow  temperature  and  chemical  composition  were  varied  to  change  the  mutual 
effect  between  mixing  and  chemical  time  scales,  and  hence  the  Damkohler  number.  The  inflow 
chemical  variation  and  the  Damkohler  number  were  found  to  affect  both  the  flame  stabilization 
and  the  lift-off  characteristics.  The  stagnation  point  position  in  the  wake  of  the  bluff-body 
and  the  lift-off  height  were  monitored  for  all  three  test  cases  to  assess  the  Damkohler  number 
similarity. 

3.1.2  Numerical  Formulation 

Filtered  Flow  Equations  -  All  simulations  reported  here  were  performed  using  the  Stan¬ 
ford  high-order  finite  difference  NGA  solver  (Desjardins  et  al.  (2008)).  The  numerical  schemes 
used  for  the  spatial  discretization  were  described  in  Desjardins  et  al.  (2008)  and  were  based 
on  the  work  of  Morinishi  et  al.  (1998).  The  describing  equations  were  in  the  low  Mach  num¬ 
ber  formulation.  The  scheme  was  staggered  in  space  and  time  (Pierce  &  Moin  (2001)).  The 
staggering  allowed  for  secondary  conservation,  which  guaranteed  the  suppression  of  numerically 
caused  instability.  Hence,  no  artificial  damping  was  required,  which  allowed  for  high  accuracy 
large-eddy  simulations.  In  this  project,  second-order  time  and  space  discretization  were  used. 
A  semi-implicit  iterative  technique,  described  below,  was  employed  to  integrate  the  system  of 
equations  using  the  fractional  step  method.  More  details  about  the  numerical  methods  can  be 
found  elsewhere  (Pierce  &  Moin  (2001),  Desjardins  et  al.  (2008)). 

By  applying  a  spatial  filter,  denoted  here  by  an  overline,  to  the  Navier  Stokes  equations  and 
using  the  Favre  decomposition,  defined  for  a  general  variable  (f)  as  (f)  =  the  continuity  and 
the  momentum  equations  can  be  written  as: 
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(1) 


dp  dpuj-  = 
dt  dxj  ’ 

dpui  dpuiUj  dp  duij  drjj 

8t.  dxj  dxi  dxj  dxj 

Here,  p  is  the  density,  U{  is  the  velocity,  p  is  the  pressure,  t  is  the  time,  Xi  are  the  spatial 
coordinates,  and  the  filtered  shear  stress  is  written  as 


aij 


2p  ( Sij 


(2) 


where  Sij  =  \  [jpp:  +  pjpj  is  the  strain  rate  tensor  and  5ij  is  the  Kronecker  symbol.  The 
summation  over  repeated  indices  is  implied,  and  p  is  the  dynamic  viscosity.  In  the  above 
equations,  the  effect  of  the  subgrid  scales  on  the  resolved  scales  is  represented  by  the  sub¬ 
filter  stresses  Tij  =  p  ( upuj  —  UiUj).  The  sub-filter  stresses  are  modeled  using  an  eddy  viscosity 
assumption  (Moin  et  al.  (1991))  as 


Uj  —  2pt 


(3) 


The  sub-filter  eddy  viscosity  pt  and  Tkk  are  modeled  as 


pt  =  ClxpA2\S\ , 

rkk  =  2CKpA2\S\2 .  (4) 

In  these  equations  A  is  the  filter  width,  |£|  =  \J 2 SijSij,  and  C ^  and  Cj<  are  computed  dynamically 
based  on  the  approach  of  Moin  et  al.  (1991).  The  density  p  is  calculated  from  the  combustion 
model  described  in  the  next  section. 

Flamelet  Progress  Variable  Approach  -  The  flamelet/progress  variable  (FPV)  approach 
developed  by  Pierce  &  Moin  (2004)  is  based  on  the  flamelet  concept  (Peters  (1984)).  This 
approach  assumes  that  the  chemical  time  scales  are  short  enough  that  chemical  reactions  occur 
in  thin  layers  around  stoichiometric  conditions  (Pitsch  (2006)).  The  flamelet  approach  relates 
the  species  mass  fractions  and  the  energy  to  the  mixture  fraction  through  the  flamelet  equations, 
which  are  given  as 

-  P2dF  =  U  •  (5) 

where  Z  is  the  mixture  fraction,  =  (T,  yi)T  is  a  vector  that  contains  the  species  mass  fractions 
Hi  and  the  temperature  T.  Co  is  the  vector  of  chemical  source  terms  for  the  scalars  T,  and  x  is 
the  scalar  dissipation  rate. 

In  the  FPV  model,  solutions  of  Eq.  5  are  tabulated  and  parameterized  by  the  mixture  fraction 
and  a  reactive  scalar  A.  As  a  result,  a  tabulated  form  T  of  the  solutions  of  Eq.  5  can  be  written 
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as  $  =  (Z,  A).  The  parameter  A  is  defined  through  the  progress  variable  C,  which  is  defined 

here  as  the  sum  of  the  mass  fractions  of  CO2,  CO,  H2O,  and  H2.  The  definition  of  A  is  discussed 
in  more  detail  in  Ihme  et  al.  (2005),  but  the  most  important  feature  is  that  it  is  defined  to  be 
independent  of  mixture  fraction.  The  mean  quantities  then  can  be  computed  from 

noo 

$  (Z,  A)  P  (Z,  dZdX ,  (6) 

where  P  (Z,  is  the  joint  probability  density  function  (pdf)  of  Z  and  A.  Since  A  and  Z  are 

defined  to  be  independent,  their  joint  PDF  can  be  obtained  from  the  product  of  the  marginal 
PDFs. 

The  mixture  fraction  marginal  PDF  P  (Z)  is  assumed  to  be  a  beta  distribution  (Pierce  & 
Moin  (2004))  and  the  reactive  scalar  PDF  a  delta  function  at  its  mean.  These  assumptions  lead 
to  a  representation  of  the  flamelet  library  (J7)  as 

$  =  Pl(z,Z^,  a).  (7) 

Writing  Eq.  7  for  the  progress  variable  as 

C  =  TX~(ZX\  A),  (8) 

the  reaction  progress  parameter  A  can  be  replaced  by  the  progress  variable.  Hence,  using  Eqs.  7 
and  8,  the  flamelet  library  can  be  written  as 

$  =  (9) 

Then,  three  parameters  have  to  be  computed,  namely  the  filtered  mixture  fraction  Z .  the  sub¬ 

filter  scalar  variance  Z "2,  and  the  filtered  progress  variable  C.  Here  the  scalar  variance  Z"2  is 
evaluated  using  an  algebraic  model  (Pierce  &  Moin  (1998)).  The  model  assumes  homogeneity 
and  local  equilibrium  for  the  sub-grid  scales  and  is  given  as 

p2r2  =  CzA2p\VZ\2,  (10) 

where  the  mixture  fraction  coefficient  Cz  is  computed  dynamically.  The  second  parameter  is 
the  filtered  mixture  fraction  Z,  for  which  a  transport  equation  is  solved 

^  +  V-  (puZ)  =V-  (p{Dz  +  DT)VZ)  .  (11) 

The  mixture  fraction  diffusivity  Dz  is  obtained  from  Eq.  (6),  and  Dt  is  the  turbulent  eddy 
diffusivity  given  by  Pierce  Sz  Moin  (1998) 


Dt  =  CtA2\S\  ,  (12) 

where  Ct  is  computed  in  a  dynamic  procedure  (Moin  et  al.  (1991)).  Finally,  the  filtered  progress 
variable  is  obtained  by  solving  the  transport  equation  of  C  as: 
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+  v  •  (puC)  =  V  •  (p(Dc  +  Dt )  VC)  +  Wc  ,  (13) 

where  Dq  is  the  progress  variable  diffusivity,  Wc  =  is  the  Favre-averaged  net  mass 

production  rate  of  C,  and  me  is  the  net  production  rate  in  (kg/m3s). 

The  flamelet /progress  variable  model  has  been  validated  in  several  a  priori  and  a  posteriori 
studies ,  including  applications  to  predict  NOx  (Ihme  &  Pitsch  (2008a))  and  local  extinction 
and  reignition  (Ihme  et  al.  (2005);  Ihme  &  Pitsch  (20086, c)).  Next  we  reformulate  the  modeled 
governing  equations  in  the  non-dimensional  form  to  derive  similarity  conditions  and  the  role  of 
the  Damkohler  number. 


Non-Dimensional  Governing  Equations  -  In  the  following,  we  considered  a  flame  stabilized 
on  a  bluff-body  flame  holder.  The  oxidizer  stream  was  vitiated,  and  the  fuel  was  injected  through 
the  flame  holder  in  cross-flow.  To  derive  the  non-dimensional  governing  equations,  the  following 
relations  were  defined,  where  the  subscript  (oo)  referred  to  the  reference  values  at  the  vitiated 
inflow  conditions  and  the  superscript  (*)  referred  to  the  non-dimensional  values: 


The  reference  length  scale  was  chosen  to  be  a  characteristic  length  of  the  flame  holder  Db-  The 
scalar  dissipation  rate  x  was  non-dimensionalized  by  its  value  at  the  quenching  conditions  \q- 
The  mixture  fraction  as  a  normalized  quantity  needed  no  further  normalization.  The  progress 
variable  mass  fraction  in  the  vitiated  coflow  was  chosen  as  a  reference  value  for  the  reaction 
progress  variable,  CTe f  =  Goo-  Another  reasonable  choice  was  to  use  the  progress  variable  value 
at  complete  conversion  of  fuel  and  oxidizer  to  reaction  products.  That  definition  could  be  used 
if  the  more  general  case  of  non-vitiated  flow  was  considered. 

By  substitution  of  Eqs.  (14)  into  the  modeled  form  of  Eqs.  (1),  the  following  non-dimensional 
form  was  obtained: 


gg  ggg  n 

dt*  dx*j 

dp*d*  dp*u*u*  =  _dp *_  _  dr 

dt*  dx*  dx*  Reoo  dx*  dx* 

J  L  J  J 


(15) 


Here,  the  Reynolds  number  of  the  vitiated  inflow  was  given  by  Re^  =  PooU°°Db .  Following  the 

M  oo 

same  procedure  with  Eqs.  (11)  and  (13),  the  filtered  progress  variable  and  the  filtered  mixture 
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fraction  equations  in  the  non-dimensional  form  became 


dp*Z 

dt* 


+L(p’~<z 


1  d 

Reoo  Scoo  dx* 


(16) 


dp*C* 

dt* 


+ 


d 

dx* 


— 

P  «i 


C‘ 


1  d  (_  ,dc*\  _d_(-D*dC^\  WcDb 

Reoo  Scoo  dx*  ^  c  dx*  J  +  dx*  ^  T  dx*  J  +  PooRfR  ’ 

(17) 


where  Scoo  =  — tw —  was  the  Schmidt  number  of  the  vitiated  inflow.  Further  treatment  for  the 
progress  variable  source  term  will  be  discussed  in  the  next  subsection.  Finally,  the  state  rela¬ 
tionship  for  the  thermodynamic  and  chemical  variables  was  computed  from  the  non-dimensional 
equation  of  state: 

**  =  rz(z,c*)  .  (is) 

Irrespective  of  the  specific  values  of  the  parameters  appearing  in  the  problem,  Eqs.  (15)-(18) 
would  have  the  same  solution  if  all  non-dimensional  groups  appearing  in  the  non-dimensional 
equations  had  the  same  value  and  if  the  non-dimensional  boundary  conditions  were  the  same. 
The  inflow  boundary  conditions  for  the  cross  flow  and  the  jet  were  uniform  with  bulk  velocities 
Uoo  and  Uj,  respectively,  where  subscript  J  stood  for  conditions  in  the  fuel  jets;  therefore,  the 
inflow  boundary  conditions  for  the  vitiated  cross  flow  in  the  non-dimensional  form  were: 


Roc)  =  1-0 
Roc}  =  o.o 
Roc}  =  o.o 
R  =  i.o 

Z  =  0.0. 


Similarly,  the  non-dimensional  jet  inflow  conditions  were: 


Rj)  =  o-o 
Rj>  =  o.o 


C*  =  0.0 
Z  =  1.0. 


(19) 


(20) 


where  r 


was  the  momentum  ratio. 


The  above  equations  showed  that  the  non- 


dimensional  solution  p* ,  u*,  p*,  Z,  C*  at  t* ,  x*  was  described  by  the  non-dimensional  groups 
Reoo,  Scqo,  r,  and  one  more  that  came  from  the  chemical  source  term  in  Eq.  17.  Assuming  that 
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the  Reynolds  number  was  large,  that  molecular  Schmidt  number  effects  were  unimportant,  and 
keeping  r  fixed,  the  only  free  parameter  was  the  Damkohler  number1,  which  was  discussed  in 
the  following  section. 

Damkohler  Number  Similarity  -  The  characteristic  Damkohler  number  was  defined  as  the 
ratio  of  the  characteristic  flow  or  mixing  time  scale  to  the  characteristic  chemical  time  scale 

Da  -  Tflow’ref  ,  (21) 

7chem,ref 

where  the  flow  time  scale  was  defined  as  tqow  =  S32-,  here  based  on  the  bluff-body  width  D b  and 
the  cross  flow  velocity  Uoq.  The  chemical  time  scale  could  be  defined  in  many  different  ways.  For 
single  step  chemistry,  it  typically  was  defined  as  the  chemical  source  term  divided  by  a  reactant 
concentration.  In  previous  work,  Pitsch  &  Fedotov  (2001)  and  Pitsch  et  al.  (2003)  derived  an 
expression  for  the  non-dimensional  source  term  for  single-step  chemistry.  They  showed  that  the 
source  term  and  equivalently  the  chemical  time  scale  were  functions  of  unique  non-dimensional 
quantities  that  depended  on  the  frequency  factor  of  the  one-step  mechanism  and  the  chemical 
and  thermal  boundary  conditions;  therefore,  mathematical  similarity  could  be  achieved  for  such 
a  case.  In  the  interaction  of  realistic  chemistry  with  flow  and  transport  processes,  different 
combustion  phenomena  occurred,  such  as  auto-ignition,  premixed  flame  propagation,  and  non- 
premixed  combustion.  Because  different  chemical  reactions  might  be  important  for  each  of  these 
processes,  a  unique  definition  of  the  relevant  chemical  time  scale  associated  with  a  particular 
reaction  was  not  possible.  Each  chemical  reaction  had  its  own  time  scale  and  Damkohler  num¬ 
ber.  As  an  example,  if  mathematical  similarity  was  to  be  achieved  at  a  higher  flow  velocity, 
which  implied  a  smaller  flow  time  scale,  the  chemical  time  scales  needed  to  be  reduced  such  that 
the  same  overall  Damkohler  number  was  achieved.  Typical  ways  to  change  the  chemical  time 
scales  were  changes  in  the  temperature,  the  equivalence  ratio,  or  dilution.  It  was  obvious  that 
not  all  chemical  time  scales  could  be  changed  in  the  same  way,  because  the  chemical  time  scales 
depended  on  both  temperature  through  activation  energy  and  species  concentrations,  the  vari¬ 
ation  was  different  for  each  reaction;  therefore,  only  one  of  many  chemical  time  scales  could  be 
adjusted.  The  choice  of  this  time  scale  should  be  related  to  the  physics  of  the  dominating  com¬ 
bustion  phenomenon.  For  example,  in  premixed  flames,  the  flame  propagation  mechanism  was 
important,  while  for  auto-ignition  stabilized  flames,  initiation  reactions  might  be  rate-limiting. 

Figure  25  shows  the  steady  solutions  of  the  flamelet  equations  for  varying  vitiated  air  temper¬ 
atures  and  stoichiometric  scalar  dissipation  rates.  As  the  vitiated  inflow  temperature  increased, 
radical-producing  reactions  were  accelerated,  and  chemical  time  scales  generally  decreased,  al¬ 
lowing  for  a  higher  quenching  scalar  dissipation  rate.  The  effect  on  scalar  dissipation  rate  was 
shown  in  Figs.  26  and  27.  Figure  26  shows  the  normalized  chemical  source  term  of  the  progress 

xNote  that  a  low  Mach  number  scaling  was  used  for  the  pressure,  which  implied  that  the  Mach  number  was 
low  enough  that  compressible  effects  were  unimportant.  This  was  the  case  for  the  conditions  considered  here; 
however,  the  assumption  that  the  effect  of  compressibility  could  be  neglected  in  the  assessment  of  static  flame 
stability  in  augmentors  needed  future  validation 
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Stoichiometric  Scalar  Dissipation  Rate,  x  [  1  /s] 


Figure  25:  Effect  of  the  vitiated  flow  temperature  (T00)on  quenching  and  ignition  characteristics 


variable  as  a  function  of  the  non-dimensional  scalar  dissipation  rate  for  the  same  variation  in 
vitiated  inflow  temperatures  as  in  Fig.  25.  The  inverse  of  a  chemical  time  scale  defined  as 


7~chem 


PooCref 


WC,  st 


(22) 


is  shown  in  Fig.  27.  As  the  dissipation  rate  increased,  which  corresponded  to  a  decrease  in  the 
flow  time  scale,  the  chemical  source  term  went  up,  and  the  chemical  time  scale  decreased.  The 
reason  was  that,  while  the  stoichiometric  temperature  decreased,  which  lead  to  an  increased 
chemical  time  scale,  stronger  mixing  increased  the  stoichiometric  concentrations  of  the  fuel, 
oxidizer,  and  related  products,  which  made  reactions  faster.  The  resulting  effect  was  a  net 
decrease  in  the  chemical  time  scale,  until  the  effect  of  decreasing  temperature  started  to  dominate 
the  chemical  source  terms,  which  ultimately  lead  to  extinction.  The  balance  between  chemistry 
and  transport  at  extinction  was  particularly  important,  since  it  provided  the  largest  possible 
chemical  source  term  and  the  smallest  possible  chemical  time  scale.  Because  the  extinction 
process  was  important  for  static  stability,  the  minimum  chemical  time  scale  corresponding  to 
the  maximum  chemical  source  term  should  be  used  here  to  define  the  Damkohler  number  of  the 
reaction  progress  variable  Dac.  With 


and  with 


Ahem,  re  I' 


Poo  Cref 
Wc,  max 


W£  = 


WC 

Wc,  max 


(23) 


(24) 


48 


1000 


Stoichiometric  Scalar  Dissipation  Rate,  xst  [1  /s] 


Figure  26:  Variation  of  normalized  chemical  production  rate  of  progress  variable  with  scalar 
dissipation  rate 


(b)  Projection  into  Xst  -  TchJm  plane 


Figure  27: 


Variation  of  chemical  time  scale  along  S-shaped  curve 
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Figure  28:  Computational  mesh  for  LES  of  the  flame  holder  geometry 


Eq.  (17)  can  be  rewritten  as 


dpC 

~df 


1  d 

Ee^o  Scoo  dxi 


+  W*cDac , 


where  Dac  was  defined  according  to  Eq.  (21). 


(25) 


3.1.3  Numerical  Method 

In  the  numerical  simulations,  the  bluff-body  geometry  was  represented  using  an  immersed  bound¬ 
ary  (IB)  technique,  as  described  by  Kang  et  al.  (2009).  The  IB  method  was  based  on  the  re¬ 
construction  of  the  velocity  components  on  the  faces  of  the  staggered  interface  cell  by  linear 
interpolation  from  the  adjacent  fluid  cells  solved  by  the  Navier-Stokes  solver.  These  recon¬ 
structed  velocities  then  were  used  as  boundary  condition  for  the  solver  around  the  IB. 

The  bluff-body  geometry  and  the  grid  are  shown  in  Fig.  28.  The  computational  grid  consisted 
of  110  X  40  x  76  mesh  points  in  the  flow  direction,  in  the  lateral  direction,  and  across  the  bluff 
body.  The  domain  was  assumed  periodic  in  the  lateral  direction,  and  the  mesh  was  refined 
in  the  regions  around  the  bluff-body.  Methane  jets  at  300  K  were  injected  from  two  circular 
ports  located  0.5  inch  upstream  the  base  of  the  bluff-body.  A  flow  of  vitiated  hot  mixture 
was  injected  in  the  stream-wise  direction  at  an  equivalence  ratio  of  =  0.67.  The  jet  and 
the  vitiated  mixture  mass  flow  rates  were  varied  to  change  the  effective  flow  time  scales,  while 
keeping  the  momentum  ratio  the  same. 

To  assess  and  validate  the  proposed  method  to  achieve  Damkohler  number  similarity,  three 
different  test  cases  were  chosen.  The  flow  conditions  for  the  three  test  cases  were  summarized 
in  Table  1. 
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Uoo 

Uj 

TociK) 

^OO 

Tflow,  ref 

T"chem,ref 

Dac 

RGqq 

SCoo 

StRZ 

Case  1 

3.6 

2.23 

1100 

1.4  •  10-4 

0.0071 

0.0011 

6.429 

653 

0.16 

0.20 

Case  2 

7.8 

4.43 

1300 

1.8  •  10-4 

0.00327 

0.00051 

6.429 

1101 

0.16 

0.18 

Case  3 

22.7 

9.75 

1700 

2.6  •  10“4 

0.00118 

1.74E-4 

6.429 

2218 

0.16 

0.22 

Table  1:  Flow  test  conditions,  where  Too  was  the  temperature  of  the  vitiated  hot  stream  [K], 
Voo  ws  the  inflow  kinematic  viscosity  [m2/s],  Too  was  the  vitiated  air  bulk  velocity  [m/s],  Uj  was 
the  fuel  jet  velocity  [nr/s],  Tflow  was  the  characteristic  flow  time  [s],  rchem  was  the  characteristic 
chemical  time  [s],  Dac  was  the  characteristic  Damkohler  number,  Reoo  and  Scqo  were  the  vitiated 
hot  flow  Reynolds  number  and  Schmidt  number,  and  StR,z  =  fRzDoo/Uoo  was  the  Strouhal 
number  of  the  oscillation  of  the  recirculation  zone,  where  /rz  was  the  corresponding  frequency. 


The  conditions  of  the  first  test  case  matched  the  experimental  setup,  while  for  the  other 
two  cases,  the  flow  and  the  chemical  conditions  were  varied  such  that  the  Damkohler  number 
was  the  same.  The  first  test  case  represented  low-speed/slow  chemistry,  while  the  second  and 
the  third  case  represented  high-speed/fast  chemistry.  The  inflow  composition  was  computed 
as  the  equilibrium  composition  of  the  methane-air  mixture  at  the  inflow  equivalence  ratio  and 
temperature.  These  values  then  were  used  to  set  the  boundary  conditions  for  the  flamelet 
equations. 

3.1.4  Results 

The  instantaneous  recirculation  zone  (RZ)  and  flame  structure  are  shown  in  Fig.  29.  The  three 
flames  showed  similar  flame  and  RZ  structure.  The  recirculation  zone  was  found  to  extend  over 
the  bluff-body  base  and  was  surrounded  by  the  flame  surface.  The  effect  of  the  outer  shear  layer 
was  shown  to  wrinkle  the  flame  surface,  causing  shedding  of  mixed  burned/unburned  gases. 

Figure  30  shows  the  mean  non-dimensional  axial  velocity  distribution  at  different  heights 
above  the  bluff-body  base.  Downstream  of  the  bluff  body,  a  recirculation  zone  was  established. 
The  RZ  extended  to  a  distance  of  around  50  mm  above  the  base  and  around  20  mm  from  the 
centerline.  Experimental  data  were  available  at  three  locations  and  were  compared  with  the 
simulation  results.  The  PIV  results  were  ensemble-averaged  over  100  frames,  with  streamwise 
averaging  resulting  in  a  statistical  accuracy  estimated  at  1  %.  Good  agreement  was  observed 
between  the  LES  and  the  experimental  data;  however,  the  measured  velocity  profiles  were  slightly 
asymmetric.  This  asymmetry  in  the  experimental  results  was  due  to  the  fact  that  the  inflow 
field  was  not  perfectly  uniform.  The  solutions  of  all  three  computed  cases  collapsed  onto  each 
other,  which  confirmed  the  Damkohler  number  similarity  of  the  three  flames. 

A  comparison  of  the  normalized  mean  temperature  is  shown  in  Fig.  31.  The  three  test  cases 
exhibited  the  same  flame  structure;  however,  at  the  first  location,  the  test  case  3  flame  showed 
a  small  deviation  at  the  edge  of  the  shear  layer.  The  mean  flame  and  RZ  structure  showed  a 
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(a)  Case  1 


(b)  Case  2 


(c)  Case  3 


Figure  29:  Iso-surfaces  of  negative  axial  velocities  (blue)  combined  with  temperature  iso-surfaces, 
T  =  1400  K  for  case  1,  T  =  1450  K  and  1550  K  for  case  2,  and  T  =  1850  K  and  1950  K  for  case 
3,  all  colored  by  progress  variable. 

mathematically  similar  solution,  again  underlining  the  applicability  of  the  proposed  Damkohler 
number  similarity  for  the  given  conditions  of  the  three  studied  flames. 

As  the  studied  flames  showed  similar  RZ  structure,  they  also  showed  similar  flame  stability 
features.  The  flame  stabilization  point  was  defined  as  the  closest  point  to  the  bluff-body  where 
reaction  happens.  Since  this  flame  was  wake  stabilized,  the  wake  structure  controled  the  flame 
stability  characteristics.  Furthermore,  as  long  as  the  chemical  time  scale  was  short  enough  to 
compete  with  residence  time,  the  flame  would  stabilize  and  would  not  blow  off.  The  shorter  the 
chemical  time  scale,  the  higher  the  flow  speed  was  needed  to  stabilize  the  flame  at  the  same 
height.  To  further  assess  the  validity  of  the  current  approach,  the  flame  height  above  the  base 
was  measured  and  compared  with  experimental  data.  By  analogy  with  the  experiment,  the 
flame  height  was  defined  as  the  first  point  above  the  base  (at  the  median  plane),  where  burning 
stoichiometric  conditions  were  found.  The  flame  height  was  computed  as  the  average  from  many 
snapshots  obtained  over  a  long  compute  time.  The  averaged  computed  values  were  around  10.5 
mm,  11.6  mm,  and  13.2  nun  for  the  three  test  cases,  respectively,  which  were  comparable  to  the 
12  mm  reported  in  the  experiment. 

Next  the  flame/wake  dynamics  were  studied  to  illustrate  the  similarity  in  the  overall  behavior. 
In  the  low-speed/slow-chemistry  flame,  case  1,  the  RZ  tip  was  found  to  oscillate  at  a  frequency  of 
28  Hz,  which  was  induced  by  the  vortex  shedding.  Figure  32  shows  eight  instantaneous  snapshots 
at  the  plane  Y  =  10  mm  for  case  1  over  a  complete  cycle.  The  temperature  distribution  was 
combined  with  the  contour  of  zero  axial  velocity  to  show  the  extent  of  the  recirculation  zone. 
The  flame  boundaries  were  shown  by  plotting  the  contour  of  the  stoichiometric  mixture  fraction 
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Figure  30:  Time- averaged  non-dimensional  axial  velocity  /U0 0  at  different  heights  X  above 

the  bluff-body  base.  Experimental  PIV  data  for  test  case  1  are  shown  as  •,  simulations  results 
are  (— )  for  test  case  1,  (★)  for  test  case  2,  and  (o)  for  test  case  3. 
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Figure  31:  Time  averaged  non-dimensional  temperature  0  =  ({T  J  —  T00)/(Tst  —  Too)  at  different 
heights  above  the  bluff  body  base,  (— )  for  test  case  1,  (*)  for  test  case  2,  and  (o)  for  test  case 
3.  Ts t  is  the  maximum  burning  temperature  at  stoichiometric  conditions 
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(Z  =  0.021).  As  the  RZ  tip  oscillated,  the  strain  rate  changed  across  the  flame  surface,  which 
was  pushed  outwards  into  the  high-speed  regions  of  the  flow.  Temperature  variations  showed 
up  at  the  edge  of  the  RZ,  indicating  high  dissipation  rates  at  these  locations.  Because  of  the 
shedding,  alternative  regions  of  hot  and  fresh  gases  were  separated  at  the  end  of  the  RZ  and 
convect  downstream.  These  regions  could  grow  in  size  as  they  progress  downstream  or  they 
could  remix  with  hot  gases  and  burned  out.  The  RZ  provided  sufficient  residence  time  to  mix 
and  combust  the  fresh  gases,  which  was  shown  in  Fig.  32.  Starting  from  Fig.  32a,  a  pocket  of 
unburned  gas  appeared  on  the  right  just  outside  of  the  RZ.  Figures  32c-32e  show  how  a  part  of 
this  unburned  gas  was  entrained  into  the  RZ,  where  it  mixed  with  hot  reaction  products  and 
burned.  The  remaining  part  of  the  pocket  was  convected  downstream,  where  it  gradually  mixed 
and  burned.  While  the  RZ  oscillated  to  the  right  because  of  the  vortex  shedding,  another  pocket 
of  unburned  gases  was  entrained  into  a  region  just  left  of  the  RZ,  where  it  underwent  a  very 
similar  process  as  the  previous  pocket  on  the  right.  This  process  repeated  in  every  cycle  on  both 
sides. 

The  same  analysis  could  be  performed  for  test  case  3,  shown  in  Fig.  33.  The  flame  showed 
wake  and  flame  structure  similar  to  the  first  flame.  The  RZ  extended  around  50  mm  above 
the  bluff-body  base.  The  oscillation  frequency  was  higher  at  196  Hz,  as  the  residence  time 
was  smaller  in  this  case,  but  the  Strouhal  number  was  very  comparable  to  the  first  case.  The 
Strouhal  numbers  for  all  three  cases  were  very  close  to  0.2.  The  values  are  given  in  Table  1.  The 
table  demonstrates  that  the  overall  dynamics  of  all  three  cases  agreed  quite  well,  which  again 
confirmed  the  applicability  of  the  Damkohler  number  similarity  concept,  as  presented  in  this 
report. 

3.1.5  Conclusion 

The  existence  of  Damkohler  number  similarity  for  static  stability  in  augmentor  flows  was  demon¬ 
strated.  A  high-speed/fast-chemistry  flame  was  shown  to  have  similar  mean  velocity  and  tem¬ 
perature  fields,  lift-off  height,  and  dynamic  behavior  as  the  low-speed/slow  chemistry  flames. 
While  the  present  study  did  not  consider  actual  blow-off,  these  characteristics  have  been  used 
as  surrogates  for  static  stability.  A  characteristic  Damkohler  number  was  defined  based  on  the 
maximum  chemical  production  rate  of  the  progress  variable.  It  was  shown  that,  for  all  cases 
studied,  the  mean  velocity  and  temperature  fields,  as  well  as  the  frequency  of  the  bluff-body 
wake  oscillation,  showed  self-similar  behavior. 

3.2  Large-Eddy  Simulations  of  New  Flame  Holder  Geometries 

In  Section  2,  experimental  studies  were  used  to  assess  the  performance  of  improved  flame-holder 
designs  for  augmentors.  Here,  we  performed  numerical  simulations  of  the  same  experimental 
setup  to  assess  the  ability  of  the  models  to  predict  the  improved  stability  observed  in  the 
new  designs.  If  successful,  the  use  of  simulation  allowed  for  the  analysis  of  the  different  flow 
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Figure  32:  Instantaneous  temperature  distribution  for  case  1  combined  with  the  contours  of 
negative  axial  velocity  (red  lines)  and  the  contour  of  stoichiometric  mixture  fraction  (black  line) 
over  a  complete  shedding  cycle 
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Figure  33:  Instantaneous  temperature  distribution  for  case  3  combined  with  contours  of  negative 
axial  velocity  (red  lines)  and  the  contour  of  stoichiometric  mixture  fraction  (black  line)  over  a 
complete  shedding  cycle 
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structures  induced  by  the  different  bluff  body  shapes  and  their  impact  on  flame  stability.  The 
interaction  of  the  different  injection  holes,  for  example,  could  be  studied  in  detail  and  optimized 
using  numerical  simulations.  As  another  example,  the  presence  of  localized  regions  of  high  fuel 
concentration  that  was  expected  in  the  local  cavity  shape  designs  could  be  investigated  with 
simulations.  Thus,  simulation  could  be  a  valuable  tool  in  the  design  and  development  of  new 
bluff  body  shapes  for  improved  flame  stability. 

3.2.1  Test  Cases  and  Simulation  Setup 

In  order  to  test  the  modeling  capabilities,  we  computed  a  subset  of  the  geometries  introduced 
in  Section  2.  The  entire  three-dimensional  experimental  setup  was  modeled,  along  with  the  six 
different  injection  holes.  A  cartesian  grid  with  local  stretching  was  used  to  represent  the  setup. 
The  generated  computational  grid  is  shown  in  Fig.  34.  The  grid  size  was  192x384x192  in  the 


Figure  34:  Grid  used  for  the  computational  study.  The  reference  shape  was  modeled  in  the  case 
shown  here.  The  blue  region  represents  the  methane  fuel  emanating  from  the  injection  holes. 

directions  x,  y ,  and  z.  As  in  the  experiments,  x  was  the  direction  aligned  with  the  methane 
injection,  y  was  the  direction  of  the  vitiated  stream  flow,  and  z  was  the  third  direction.  The 
location  of  the  different  holes  can  be  seen  in  Fig.  34.  Appropriate  stretching  was  introduced  to 
resolve  the  geometric  shapes  adequately.  The  simulation  of  the  entire  channel  was  avoided  by 
trimming  the  domain  to  8  inches  wide  in  the  x  direction  and  14  inches  in  the  2  direction.  A  zero 
velocity  gradient  condition  was  applied  on  the  lateral  surfaces.  The  extent  in  the  y  direction 
was  4  inches  upstream  and  5  inches  downstream  of  the  bluff  body.  The  crossflow  had  a  bulk 
velocity  of  5  m/s  and  a  turbulence  intensity  of  8  %.  An  isotropic  turbulence  field  was  generated 
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and  injected  at  the  vitiated  air  inflow  at  the  rate  of  the  bulk  velocity.  The  temperature  of 
the  crossflow  was  set  to  1315  K  based  on  experimental  measurements.  The  composition  of  the 
crossflow  was  obtained  through  a  premixed  combustion  calculation  for  an  equivalence  ratio  of 
0.67.  The  methane  jet  had  an  initial  temperature  of  400  K  and  the  same  momentum  flux  as  the 
crossflow. 

The  experimental  cases  used  for  the  study  were  the  Reference  case,  the  D+S  (in-phase),  and 
the  T  shapes.  The  geometries  for  the  different  cases  are  depicted  in  Fig.  4.  It  was  found  in  the 
experiment  that  the  Reference  case  had  a  high  liftoff  height  of  around  37  mm.  The  D+S  shape 
had  localized  cavities,  which  caused  a  lower  liftoff  height  of  around  17  mm.  The  T  shape  also 
had  a  small  liftoff  height  of  around  22  mm.  The  cases  were  chosen  to  span  the  range  from  low 
liftoff  to  high  liftoff.  Modeling  these  geometrically  varied  shapes  was  also  a  good  test  of  the 
robustness  of  the  numerical  simulation  code. 

3.2.2  Sub-Filter  Models  and  Numerical  Schemes 

The  code  used  for  this  study  was  a  high-order  kinetic  energy  conserving  code.  The  description 
of  the  code  and  the  equations  solved  has  been  provided  previously  in  Section  3.1.2.  LES  studies 
were  performed  by  solving  the  filtered  flow  equations  that  were  presented  earlier.  A  dynamic 
Smagorinsky  model  with  Lagrangian  averaging  was  used  to  model  the  sub- filter  stresses.  The 
chemistry  in  the  simulation  was  solved  using  the  Flamelet  Progress  Variable  model  (FPV)  de¬ 
scribed  in  the  previous  section.  The  chemistry  was  parameterized  as  function  of  the  progress 
variable  (C),  mean  mixture  fraction  ( Z )  and  the  variance  of  the  mixture  fraction  ( Z //2).  The 
composition  at  zero  mixture  fraction  corresponded  to  the  conditions  of  the  vitiated  crossflow. 
A  unity  mixture  fraction  value  corresponded  to  pure  methane  fuel. 

3.2.3  Results 

The  two-dimensional  instantaneous  velocity  fields  for  the  three  different  cases  are  shown  in 
Fig.  35  for  the  Pxy_C  plane.  The  flow  features  in  the  simulation  were  similar  to  those  observed 
in  the  experiments.  A  stagnation  region  formed  upstream  of  the  bluff  body,  as  expected.  The 
methane  jets  in  the  crossflow  should  lead  to  a  stagnation  region  downstream  of  the  jets,  which 
was  observed  in  the  simulation  results.  Downstream  of  the  bluff  body,  recirculation  regions  were 
formed  for  all  three  cases.  The  recirculation  region  wes  found  to  be  closer  for  the  T  and  D+S 
shapes  than  the  Reference  case.  This  trend  also  was  seen  in  the  experiments.  Centerline  profiles 
of  vertical  velocity  downstream  of  the  bluff  body  are  shown  in  Fig.  36.  The  velocity  profiles 
compared  reasonably  well  with  the  experimental  plot  shown  in  Fig.  lOd  for  the  reference  and 
the  triangle  cases.  For  the  D+S  case,  the  velocity  was  higher  than  the  experiment,  but  a  similar 
trend  is  observed. 

The  temperature  distributions  for  the  different  cases  are  presented  in  Fig.  37.  The  stoichio¬ 
metric  mixture  fraction  iso-surface  also  is  indicated  in  these  plots  by  the  black  line.  The  mass 
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Figure  35:  Instantaneous  velocity  magnitude  field  for  the  three  different  bluff  body  shapes. 


Figure  36:  Computed  vertical  velocity  profiles  in  the  y  direction  along  the  centerline. 

fraction  of  the  CH  radical  is  shown  in  Fig.  38  to  compare  with  visual  images  of  CH  from  the 
experiments.  Flames  were  seen  for  all  cases  at  some  distance  downstream  of  the  jets;  however, 
there  were  significant  differences  from  the  experimental  results.  The  main  difference  was  that 
the  flames  form  on  the  edge  of  the  bluff  body  wake,  which  was  not  seen  in  the  experiment. 
Furthermore,  the  flame  was  attached  to  the  flame  holder  in  the  region  downstream  of  the  bluff 
body  for  all  computed  cases,  which  also  was  not  observed  in  the  experiments.  In  the  simulation, 
the  mixture  fraction  was  close  to  stoichiometric  on  the  edges  of  the  wake  regions,  which  could 
be  observed  in  Fig.  37.  Thus,  it  was  reasonable  to  see  a  high  temperature  region,  since  the  flame 
typically  had  the  highest  temperature  close  to  stoichiometric.  Such  an  attached  flame  along 
with  the  wake  edge  flame  was  seen  in  the  experiments  for  certain  design  geometries.  Examples 
were  shown  in  the  phase  I  report.  It  was  found  in  the  experiments  that  introducing  the  tur- 
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Reference  T  D+S 

Figure  37:  Instantaneous  temperature  fields  for  the  three  different  bluff  body  shapes.  The  black 
line  represents  the  iso-contour  of  stoichiometric  mixture  fraction. 
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Figure  38:  Instantaneous  CH  field  for  the  three  different  bluff  body  shapes. 
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bulence  mesh  changed  the  flame  topology  from  a  jet-stabilized  to  a  wake  stabilized  flame.  The 
effect  of  the  turbulence  grid  was  represented  by  using  isotropic  turbulence  at  the  inlet  at  the 
measured  turbulence  level;  however,  this  turbulence  grid  model  might  not  be  accurate  enough, 
and  the  turbulence  mesh  would  be  computed  explicitly  in  the  future.  It  also  was  possible  that 
the  liftoff  was  sensitive  to  other  boundary  conditions  or  modeling  assumptions  prescribed  in  the 
simulations.  We  therefore  performed  a  sensitivity  study  to  test  the  importance  of  some  of  the 
assumptions  in  the  simulations. 

Among  the  conditions  specified,  the  vitiated  flow  temperature  was  particularly  important 
for  the  flame  characteristics.  Its  effect  was  studied  here  by  varying  the  inlet  temperature  to  1200 
K  and  1100  K.  Another  assumption  made  in  the  chemical  model  used  for  the  simulations  was 
the  use  of  unity  Lewis  numbers.  The  justification  for  selecting  unity  Lewis  number  was  that  in 
turbulent  flows,  the  turbulence  was  usually  intense  enough  that  turbulent  diffusivity  dominates 
molecular  diffusivity.  Thus,  all  species  diffused  in  a  similar  way.  The  effect  of  this  assumption 
was  tested  by  performing  an  additional  simulation  using  non-unity  Lewis  numbers. 

3.2.4  Effect  of  Crossflow  Temperature 

The  crossflow  temperature  was  varied  to  1100  K  and  1200  K  to  see  the  sensitivity  of  the  flame 
location  to  this  parameter.  The  experimentally  measured  crossflow  temperature  was  1315  K. 
This  variation  was  done  for  the  reference  shape  of  the  bluff  body.  It  was  found  for  both  these 
cases  that  the  wake  edge  flames  were  still  present.  The  flame  also  was  attached  to  the  bluff 
body  for  these  cases.  Uncertainties  in  the  crossflow  temperature  therefore  could  not  explain  the 
observed  differences  with  the  experiment. 

3.2.5  Effect  of  Lewis  number  Approximation 

The  influence  of  non-uniform  diffusion  was  checked  by  performing  an  additional  simulation 
with  non- unity  Lewis  numbers.  For  this  simulation,  the  crossflow  temperature  was  set  to  the 
experimentally  measured  value  of  1315  K.  The  simulation  was  performed  for  the  reference  shape 
of  the  bluff  body.  The  results  for  the  CH  mass  fraction  and  temperature  are  shown  in  Fig.  39. 
The  Lewis  number  did  not  have  a  significant  effect  on  the  results.  Both  the  temperature  and 
CH  mass  fraction  were  similar  to  the  results  assuming  unity  Lewis  numbers. 

3.2.6  Conclusions 

The  experimental  setup  used  to  study  bluff  body  shapes  in  augmentors  was  simulated  using 
large  eddy  simulation.  The  velocity  comparisons  illustrated  that  the  simulation  had  similar 
flow  features  as  were  observed  in  the  experiments.  The  temperature  and  the  CH  mass  fractions 
showed  that  the  flame  was  attached  to  the  bluff  body  for  all  cases.  Also,  flames  formed  on  the 
edge  of  the  bluff  body  wake.  Both  these  features  were  not  seen  in  the  experiments.  The  effect 
of  boundary  condition  uncertainty  and  the  unity  Lewis  number  assumption  in  the  model  were 
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Figure  39:  Instantaneous  temperature  and  CH  field  for  the  non-unity  Lewis  number  study.  The 
black  line  represents  the  isocontour  of  stoichiometric  mixture  fraction. 


investigated  by  performing  additional  simulations  with  varied  crossflow  temperature  and  with 
the  non-unity  Lewis  number  assumption.  Both  sensitivity  studies  did  not  reveal  a  significant 
effect. 

It  was  speculated  that  the  reason  for  disagreement  was  a  higher  level  of  inflow  turbulence  in 
the  experiment  or  a  very  non-isotropic  flow  induced  by  the  turbulence  mesh.  The  cited  reason 
was  confirmed  by  several  observations:  i)  flame  stabilization  on  the  fuel  jets  was  observed  in 
the  phase  I  experiments.  After  adding  the  turbulence  grid,  the  stabilization  point  in  the  exper¬ 
iments  moved  into  the  wake  of  the  bluff  body,  which  confirmed  the  importance  of  the  incoming 
turbulence,  ii)  The  CH  chemiluminescence  plots  in  Fig.  3  showed  substantial  turbulence  in  the 
wake  region,  which  was  not  observed  in  the  simulations  and  which  probably  originated  from  the 
inflow  turbulence,  iii)  The  inflow  velocity  profile  and  Reynolds  stress  tensor  were  not  very  well 
quantified  in  the  experiment;  therefore  the  prescribed  values  had  high  uncertainty. 

In  the  future,  we  shall  use  the  AFRL  experiments  of  the  flame  holders  proposed  here  for 
further  validation. 
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